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ABSTRACT 

In this work we have carried out an in-depth analysis of the young stellar content in the W3 
CMC. The YSO population was identified and classified in the IRAC/MIPS color-magnitude space 
according to the 'Class' scheme and compared to other classifications based on intrinsic properties. 
Class O/I and II candidates were also compared to low/intermediate-mass pre-main-sequence stars 
selected through their colors and magnitudes in 2MASS. We find that a reliable color/magnitude 
selection of low-mass PMS stars in the infrared requires prior knowledge of the protostar population, 
while intermediate mass objects can be more reliably identified. By means of the MST algorithm and 
our YSO spatial distribution and age maps we investigated the YSO groups and the star formation 
history in W3. We find signatures of clustered and distributed star formation in both triggered 
and quiescent environments. The central/ western parts of the GMC are dominated by large scale 
turbulence likely powered by isolated bursts of star formation that triggered secondary star formation 
events. Star formation in the eastern high density layer also shows signs of quiescent and triggered 
stellar activity, as well as extended periods of star formation. While our findings support triggering 
as a key factor for inducing and enhancing some of the major star forming activity in the HDL 
(e.g., W3 Main/W3(0H)), we argue that some degree of quiescent or spontaneous star formation is 
required to explain the observed YSO population. Our results also support previous studies claiming 
an spontaneous origin for the isolated massive star(s) powering KR 140. 

Subject headings: ISM: clouds — ISM: individual (Westerhout 3) — infrared: stars — stars: formation 
— stars: protostars — stars: pre-main sequence 



1. INTRODUCTION 

The W3 giant molecular cloud (GMC) together with 
W4 and W5 form part of a massive star forming complex 
in the P erseus Arm. With an estimated mass of ^ 4 x 
10^ M© dMoore et al.|2007[ [Polychroni et al.|2010| and a 
size of ~ 1.5 deg^ this cloud is one of the most massive in 
the outer Galaxy. Located at a distan ce of ~ 2kpc (e.g., 
Xu et al.||2006| |Navarete et al.||2011|), W3 comprises a 



tures ranging from the largest. Galactic scales, to giant 
molecular clouds, disks, and planetary syst ems. The- 



wealth of H II regions, clusters, and active star forming 
sites with clear signatures of both low-mass and high- 
mass star formation. 

Its relatively close location and its massive star con- 
tent has made this cloud a prime target f or the study of 
cluste r and high-mass star formation (see Megeath et al. 
[20081 for a detailed review of the literature related to 
the study of W3), which contrary to the low-mass case 
is, overall, poorly understood. Massive OB (M> 8M0) 
stars are sources of UV photons, responsible for disso- 
ciating/ionizing molecules and atoms (HII regions), and 
of metals, which enrich the surrounding gas with heavier 
species. They also affect the local and global dynamical 
state of the interstellar medium (ISM) through ioniza- 
tion and radiation pressure, stellar winds, outflows and 
supernovae. These stars are therefore fundamental in 
determining the state and evolution of the ISM, as well 
as the formation, maintenance, and dissipation of struc- 
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ories such as the turbulent core model (Mc Kee fc Tan 
2003) and competitive accretion (e.g., Bonnell & Bate 
2005 ) have been proposed as possible theories explammg 
the formation of the rare and massive OB stars. How- 
ever, the question as to whether massive star formation 
is simply a scaled-up version of low-mass star formation, 
or if it is the result of a completely different process, 
remains one of the main outstanding issu es in star for- 
mation theory. Indeed, authors such as Yorke ( 1993 1 
claim that a bimodality should arise from the snorter 
Kelvin-Helmholtz timescale for massive stars, their high 
L/M ratio, and the dramatic effects their strong radia- 
tion field have on their evolution. The immediate effects 
of high-mass star activity can be observed, for instance, 
with masers, shocks, 'hot spots', and the formation of hy- 
percompact, ultracompac t and classical HII regions (e.g., 
Zinnecker & Yorke 2007 and references therein). Fur- 



thermore, while low-mass stars are known to be able to 
form in isolation, most star formation is thoug ht to oc- 
cur in clusters embedded in their parent GMCs ( Lada & 
Lada 2003). This is particularly true for massive stars, 
which have been theorize d to form exclusively in clusters 
(e.g., |de Wit et alT]|2005 l, making cluster studies crucial 
in order to understand and investigate the formation of 
massive stars. W3 is believed to contain massive stars in 
various evolutionary stages (e.g., Tieftrunk et al.||T997 ). 
The eastern high density layer (HUL) neighboring W4 
contains the most active star forming sites W3 North, 
W3 Main, W3 (OH) and AFGL 333, with signatures of 
massive star for mation in a triggered environment (e.g., 
Oey et al.|2005 ) . A rare case of spontaneous massive star 
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Fig. 1. — Greyscale Spitzer channel 1 mosaic, with labels marking the regions and key features in the W3 GMC. 



formation is also belie ved to have occurr ed in KR 140, 
west of the HDL (e.g., [Kerton et al.||2008[ ). 



1.1. Proposed Analysis: The W3 GMC 

This work is the first of a series of papers focus- 
ing on W3 and which aim to shed some light on the 
still poorly understood massive star formation process. 
With the advent of high sensitivity instruments such as 
Spitzer and the Herschel Space Observatory, this GMC 
presents an ideal opportunity to investigate the high- 
mass star/cluster star formation process and its relation 
to triggered/spontaneous star formation events, clus- 
tered, and distributed star formation. In this analysis 
we fol low up the study initiated by |Polychroni et aL] 
(2010) to characterize the pre-stellar population in W3 



by means of Spitzer data. This analysis will complement 
the Herschel data obtained under the Guaranteed Time 
Key Program HOBY^ (Herschel imaging survey of OB 
young Stellar objects), as well as an extensive molecular 
analysis (Polychroni et al., in preparation), which focuses 
on the formation, evolution, and dynamics of the largely 
unknown progenitors of massive stars and clusters in this 
cloud. 

In Section [2] we first describe the infrared data used 
to select and classify the young stellar objects (YSOs). 
Classification techniques and schemes, including a de- 
scription of the final catalog of candidate YSOs have been 
included in Section l3] The physical and observational 
properties of the YSO sample are presented in Section 

t Their clustering and spatial distribution, as well as 
e techniques employed in our analysis, are described 
in Section [5] Section |6] investigates the cluster proper- 
ties and ages across W3, and results are used to compare 
the star formation history and activity in the different 



subregions. We conclude with a summary of the star 
formation in W3. 

2. PHOTOMETRY, DATA PROCESSING & DATASETS 
2.1. Spitzer IRAC and MIPS Infrared Observations 

The main active regions of W3 were observed by 
Spitzer under two main Program Identification Number 
(PID) programs. The northern parts of the HDL com- 
prising W3 Main and W3 (OH) (Figs|l]andj2| were first 
observed in 2004 under program POoM; AFGL 333 and 
central/ western region of W3 (including KR 140 and the 
active region to its north: 'KR 140-N') were subsequently 
observed in 2007 (program P30955). A preliminary anal- 
ysis of the data ob tained during the first obse rvation has 
been presen ted in Ruch et al. (|2007), while Polychroni 
et al. (2010) combined data from both programs tor their 
analy sis ot the stellar po pulation associated with SCUBA 
cores ( Moore et al. 2007 ) . For this work we also used data 
from both programs. 

For the reduction process we downloaded the IRAC 
and MIPS 24 fj,m corrected Basic Calibrated Data 
(cBCD) from the Spitzer Archive. The IRAC cBCDs 
were produced with the S18.7.0 pipeline. P00127 and 
P30955 MIPS data were processed with the S18.13.0 and 
SIS. 12.0 pipelines, respectively. These files already in- 
clude all the pre-processing artifact mitigation correc- 
tions, including muxbleed, column pulldown/puUup and 
electronic banding. Tiles were reduced, background- 
matched and mosaicked with the MOPE>s[^ tool, which 
with the exception of the MIPS observations from pro- 
gram P30955 (KR 140) produced maps of higher qual- 
ity than the mosaicked pBCD data already provided by 
the Archive. The best available mosaics were chosen for 
our analysis. Photometry tests performed on both the 
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Fig. 2. — Greyscale Spitzer channel 4 mosaic of W3. Intensity scale has been chosen to highlight the weaker features surrounding the 
arming regions labeled in Fig[T](at the expense 
and filaments (e.g., (a, 5) = (02h 26m 57^761° 29' 45"))- 



main star forming regions labeled in Fig[T](at the expense of the latter). Details in this image include several Infrared Dark Clouds (IRDCs) 
= (02h 26m 57%7 f 



pBCDs and our own reduced mosaics shows no system- 
atic differences in channels 3 and 4, and negUgible (0.04 
mag) differences in channel 1 and channel 2, which is 
much smaller than any of our expected photometric un- 
certainties. 

2.2. Spitzer Source Extraction and Photometry 

A pre hminary list of sources was obtained using SEx- 
tractor (Bertin & Arnouts 1996 1 and a mexican hat filter, 
which was observed to perforni the best (with the high- 
est detection rate of visually confirmed sources in the 
region) in the extremely crowded regions found in W3, 
especially W3 Main and W3 (OH). We note that this 
also introduced a significant number of artifacts and false 
detections which we eliminated through a series of clean- 
ing steps, described below. SExtractor also produced 
noise and background maps, which we also checked vi- 
sually to ensure optimal source extraction. This prelim- 
inary list was subsequently fed into the point source ex- 
traction package in MOPEX (APEX) for point response 
function (PRE) fitting, which provides a more accurate 
centroid calculation as well as additional statistics for 
each source resulting from the fitting (e.g., x^, signal-to- 
noise ratio (SNR), etc). When bandmerging the IRAC 
list with the 2MASS Point Source Catalog (through the 
GATOR interface in the NASA/IPAC Infrared Science 
Archive; IRSAQ we find the accuracy of our final PRE- 
fitted coordinates to be generally better than 0".5, al- 
though in this work we allow for a more conservative 
matching radius of 2" when bandmerging catalogs at 
different wavelengths and/or instruments. The internal 
background image used in APEX for photometry was 
produced choosing the option 'SExtractor background 
method' for consistency with the previous part of this 

^ http:/ /irsa. ipac.caltech.edu/applications/Gator/ 



analysis. 

As recommended by the IRAC handbook we chose 
aperture photometry for our main photometric analysis, 
using in this case the most accurate centroids returned 
by APEX. The use of this technique also ensured consis- 
tency with the most re cent Spitzer studies of this cloud 
( [Polychroni et al.||2010[ ). 
Aperture corrections and zero points were obtained 



from the IRAC/MIPS Instrument Handbook^ Ad 



ditional corrections (e.g., pixel phase and array correc- 
tion) were applied when required. An aperture of 2 pixels 
(2". 4) with a sky annulus between 2 and 6 pixels (7". 3) 
for IRAC, and an aperture of 7".0 with a sky annulus 
between 7 — 13" for MIPS, were found to yield results 
most closely agreeing with the magnitudes provided by 
Ruch et al. J 2007 ) , who used the version of daophot 
1 Stetson|19a7[ ) modified by the Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire (glimpse). Contrary 
to standard aperture photometry, the glimpse technique 
is particularly useful for analysis in crowded fields and 
regions with variable background and we therefore 
checked the accuracy of our results by comparing the 
photometry for those YSO candidates in our list with 
counterparts in the source list provided by these au- 
thors. The root mean square (rms) difference between 
the MIPS 24/im photomet ry obtained in th is work for 



the YSO list and that from Ruch et al. (2007) was found 
to be ^ 0.6 mag. Comparison of IRAC photometry yields 
an rms difference of < 0.2 mag in all four channels, which 
is consistent with the estimated 3(7 errors and equivalent 
to the minimum signal-to- noise (S/N) required for our 
final catalog (S/N= 5; see below). 



* http:/ /ssc. spitzer. caltech.edu/irac/iracinstrumenthandbook/ 
^ http:/ /ssc. spitzer. caltech.edu/mips/mipsinstrumenthandbook 
see http://www.astro.wisc.edu/glimpse/photometry 
vl.O.pdf; http://www.astro.wisc.edu/glimpse/ 
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The same procedures applied to the IRAC long expo- 
sure mosaics were also performed on the short exposure 
images. While the final catalog is based on the long expo- 
sure maps, we used the short exposure mosaics to obtain 
replacement photometry for those sources with bad pix- 
els within their apertures or observed/suspected to be 
affected by saturation in the long exposure mosaics. 

3. STELLAR CLASSIFICATION 

3.1. General YSO Classification: Methodology and 
Techniques 

3.1.1. Protostars & Optically Thick Disks 

We aimed to provide the most reliable sample of young 
stellar objects (YSOs) in this GMC. For the main clas- 
sification in our analysis we chose the 'revis e d', up dated 
criteria in Appendix A of Gutermuth et al. (2009). The 
color and magnitude scheme developed by these authors 
includes a series of sequential steps (phases) to identify, 
clean, and classify the candidates: Phase 1) Removal of 
contaminants such as star-forming galaxies with strong 
polycyclic aromatic hydrocarbon (PAH) emission, AGN, 
unresolved knots of shock emission, and PAH-cmission- 
contaminated apertures. This step also includes the 
first separation of YSOs by means of the four IRAC 
bands. Phase 2) A search for additional YSOs based 
on 2MASS photometrjj^ Phase 3) Identification and 
re-classification of previously identified sources with suit- 
able MIPS 24 /im photometry. When re-classifying pho- 
tospheric sources into 'transition disk' objects (those 
with significant 24 /im emission; included within the 
Class II category) we required sources to have been classi- 
fied as photospheric in both previous phases. To dered- 
den the magnitudes we u sed the extinction maps an d 
methodology described in Rowles & Froebrich (2009). 
The visual extinction map was transfor med to a median 
Ah map using the extinction law from Mathis (1990). 
We changed extinction in the 2MASS bands to e xtinction 
in the IRAC channels using the numbers from |Flaherty| 
eraL] ( p007| . 

i'he above classification was complemented and cross- 
checked with: i) the 'red source' classification scheme 



phase from Robitaille et al. 

6 



from Robitaille et al. (2008), which should include all 
Class O/I and sev eral Class II so u rces; and ii) the 

'stage 

Monv/M^ > 10 " yr 
and Mdisk/M^ > 10" 
and Mdisk/M,^ < 



( |2006| (Stage O/I: 
6yr- 

6 ,,^-1 



JM, < 10- 



Stage II: M< 
6; Stage III: M.^./M, < 10 "yr 
10~^). We used this last scheme to 
compare the above observational classification with an 
alternative method based on intrinsic physical proper- 
ties (e.g., mass accretion rate and disk mass). This last 
analysis was carried out by studying the position of each 
YSO in the color-color diagrams (CCDs) with respect to 
the limits marking the areas where most of the sources 
of one part icular 'stage' are predicted to fall ( (Robitaille 
eFaL||2006l ). 



3.1.2. Pre-Main Sequence Population With Optically Thin 

Disks 

O ne o f the equations of P hase 2 ( |Gutermuth| 
^ . , , „„„„ should read: 



Separation and classification of pre-main sequence 
stars (PMS) with optically-thin disks is a particularly 
complicated process due to their similarity (in infrared 
color) with more evolved reddened main sequence and 
giant stars (photospheric-dominated). These transition 
objects are however essential to fully understand the dif- 
ferent stages in star formation. In an attempt to estimate 
the population of sources with weak infrared excess and 
other PMS stars that may have been missed with the 
above color classification, we first excluded those objects 
in our source list already classified as Class O/I and II. 
We then used the 2MASS catalog and a process simi- 



lar to that used in Kerton et al. ( 2008 1 , who attempted 



to separate the YSO population guided by a sample of 
known low-mass T-Tauri and intermediate-mass Herbig 
Ae/Be (HAeBe) stars. In this work (see Section |4) we 
show that this method can only be applied successfully 
once the 'younger' YSOs have been identified using ad- 
ditional data (e.g., IRAC). 

To investigate the possibility of missed candidates 
from t he PMS population we firs t chose a sample of T- 
Tauri (Kenyon fc Hartmann' 1995 1 and HA eBe star s with 
know n distances (I'lnkcn.zgljgr fc Mundt|1984, The et al. 
1994| Mendigutia et al. |2Ullj and references therein). 
All sources where checked with SiMBAEp^ keeping vari- 
able, emission and pre-main sequence stars and reject- 
ing those classified as double/multiple systems, low mass 
stars and brown dwarfs. Infrared photometry was ob- 
tained by matching the sample with the 2MASS Point 
Source Catalog through the GATOR interface. Infrared 
photometric systems were converted to the 2MASS sys - 
tem using the transformations from Carpenter (2001). 
All magnitudes were shifted to a distance of 2 kpc for W3 
with the inclusion of interstellar e xtinction by m eans of 
the Ay— distance relation from |Indebetouw et al. ( 2005) . 
Figures [s] and |4] show the color-color diagram (CCD) 
and color-magnitude diagram (CMD) for the T-Tauri 
and HAcBc samples shifted to the distance of W 3 . The 
CCD shows the T-Tauri locus from [Meyer et al.|(|1997l 



and the main sequence and giant branch from 
(jl^983) including interstellar reddening. Reddening vec 
tors for an Ay — 10 have been included for an 06 V 
M8V, M2V, and an M6III star. The CMD shows so- 



lar metallic ity isochrones (Marigo et al. 
et al.||2010[ ) at log(tyr-i) 
tance! ^--^--J j-^^- 

for an 



2008 Girardi 



7, 8, and 9 for the same dis- 
TTie dashed-dotted line is the reddening vector 
AO star with Ay = 10 at d= 2 kpc, applied us- 
ing the redden ing law and the Ay/Aj conversion from 
Mathis (1990). The values used for extinction conver- 
sions between 2MASS bands are consistent with those 
from 'Indebetouw et al. ("2005) up to one decimal place 



|et al.| |2008| |Gutermuth et al.l |2009^ 

E[3.6]-[4.5]/E_ff-Jf = (Ejif-if /Ejf _[4.5])~-^-(Ejj_j; 



taking into account uncertainties, which is a negligible 
difference compared to the expected uncertainty in the 
transformation from optical to infrared extinction. 

Figure [3] shows T-Tauri stars lying mainly above the T- 
Tauri locus. Many are withi n the redde ning band formed 
by the reddening vectors (Mathis 1990) of an 06 and 
~M2 main sequence stars (with a tail extending into the 
HAeBe region, to the right of the reddening band, and 
following the direction of the T-Tauri locus). The wide 
distribution implies variable amounts of extinctions to- 



http : // simbad . u-strasbg.fr / simbad / 
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TABLE 1 

Selection Criteria for, T-Tauri and HAeBe stars for 



Fig. 3. — Color-Color Diagram showing the T-Tauri and HAeBe 
samples shifted to d= 2kpc. Solid lines mark the main sequence 
and giant branch from |Koornneef| l |1983[ l , also shifted to a distance 
of 2 kpc. Dash-dotted lines are reddening vectors for an additional 
Av = 10 for an 6V. MSV. M 2V. and M6III star using the 
extinction law from IMathis (1990). Dashed line marks the locus 
of T-Tauri stars from Meyer ct al. (1997) at the same distance. 
Vertical solid line marks the bluest [H-KJ color accepted for PMS 
classification. 




TTouri O 
Herbig ^ 



Fig. 4. — Color-Magnitude Diagram for the T-Tauri and HAeBe 
amples at d= 2 kpc Solid lines are sol ar metallicity isochrones 
Marigo et ar]|2008| [Cirardi et al.||2010| log(t yr~i)=7, 8 and 9. 
Jasned-dotted line is tne reddening vector for an ~A0 star at this 
distance with Av = 10. Horizontal solid line marks the magnitude 
limit separating T-Tauri and HAeBe candidates. Vertical solid line 
like in Fig[3] 

ward these sources and variable disk emission. A large 
proportion of T-Tauri stars are indistinguishable from 
main sequence stars with just interstellar redde ning or 
are consist ent with weak-emission T-Tauri stars jMeyer 
et al.|1997|). The maximum Ay in the maps from'Kowles 
fc I'Voeb rich] (]2009) is about - 9.5 - 10 in the region con- 
taining W3 Main/(OH), and ~ 7 for those comprising 
KR 140 and AFGL 333. Thus there wiU be considerable 
extinction of objects in the W3 field, and so color iden- 
tification of this type of T-Tauri star without the aid of 
spectroscopic data will be severely contaminated. 

In consequence, we chose a more conservative approach 
and selected our T-Tauri sample by requiring these PMS 
stars to be reddened enough to lie above the T-Tauri 
locus, with colors sati sfying H-K> 0.7 (similar to that 
used in |Kerton et al.| (2008) for KR 140). The above 



D= 2 KPC 


T-Tauri 




[H-KJ> 0.7 






-H]> 0.59([H-K]-0.187) H 


h 0.72 




-H]< 1.55([H-K]-0.39) 4- 


0.85 




K> 12 




HAcBc 




[H-KJ> 0.7 




[J- 


H]< 1.55([H-K]-0.187) - 


0.008 




K< 12 





limits minimize the contamination from foreground or 
mildly reddened weak-emission T-Tauri stars (undistin- 



guishable from main sequence), and early type stars. 
While the color cuts could still allow for non- negligible 
contamination from late main sequence and giant stars 
with moderate reddening, the magnitude selection cri- 
terion below reduces this contamination to mainly that 
caused by highly extinguished (Ay > 7) stars of spectral 
type of ^ A or later (Figure B . 

HAeBe stars lie preferentially to the right of the red- 
dening band. All suitable candidates should therefore 
be located in this region and satisfy the condition H- 
K> 0.7. This limit aims to minimize contamination from 
reddened early type stars and luminous T-Tauri stars. To 
separate candidates in regions of the CCD populated by 
both types of PMS (i.e., outside the reddening band) we 
used the information from the CMD (Figure l4|, in which 
T-Tauri and HAeBe can be easily separated. ^11 T-Tauri 
stars have K magnitudes > 12. HAeBe stars tend to 
be brighter, approaching this limit only for 'late' stages 
reaching the main sequence, which we already discard 
with our imposed limit in the CCD to minimize contam- 
ination from reddened early type stars. The combined 
color plus magnitude condition minimizes contamination 
of the T-Tauri sample from reddened giant stars. In addi- 
tion, the color constraint imposed on HAeBe stars, which 
are mainly localized and 'isolated' outside the reddening 
band of typical stars, already minimizes the contamina- 
tion from reddened main sequence and giant stars. 

We note that this relatively simple scheme for PMS 
classification may only be applied after the Class O/I 
and Class II populations have been identified using IRAC 
data, as there is considerable overlap in the CCD and 
CMD of T-Tauri candidates with Spitzer Class II/I/O 
sources (see Section |4]). 

Our final PMS selection scheme has been summarized 
in (Table |T]). The color and magnitude selection criteria 
were applied to all sources in our initial IRAC list satis- 
fying i) the cleaning/reliability conditions in the Spitzer 
channels; ii) matched to a 2MASS source with quality 
flag better than 'D' in all 2MASS bands; and iii) classified 
as (mainly) photospheric or without a successful classi- 
ficatio n using the YSO scheme from Gutermuth et al. 
(20091 (i.e., no Class O/I, Class II, or contaminant). 

A search for additional PMS was also carried out by 
extending our analysis to 2MASS sources in the area cov- 
ered by the Spitzer survey but without a suitable IRAC 
counterpart (i.e., satisfying our initial cleaning and reli- 
ability conditions) in the short wavelength Spitzer chan- 
nels. The lack of a detection in Spitzer minimizes the 
possibility of confusion with actual embedded protostars, 
which should have been prev iously identified with th e 
color/magnitude criteria from Gutermuth et al. (2009). 
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3.2. The Spitzer Catalog 

Here we present the final products and source fists 
derived from the analysis carried out in the previous 
section. Results and statistics from our YSO detection 
and classification procedures are shown in Tables [2] and 

tFor the purpose of this paper, we define as 'YSO' 
ose Class O/I and Class II ca ndidates selected using th e 
color/magnitude scheme from Gutermuth et al. (2009). 
The sample of PMS stars, that is those additional can- 
didate young stellar objects selected using 2MASS pho- 
tometry, will include stars with optically thick disks, e.g., 
classical T Tauri stars and HAeBe stars (Class II) missed 
by the IRAC color/magnitude classification, and opti- 
cally thin disks, e.g., weak-lined T Tauri stars (Class III 
sources). An analysis of this sample will be relevant to 
investigating the 'oldest' young stellar population in W3. 



3.2.1. Statistics, Completeness & Reliability 

The YSO s earch using the color and magnitude classi- 
fication from Gutermuth et al. (2009) yielded a total of 
616, 706 (two ot which were also observed in neighboring 
AORs), and 246 YSOs in the regions surveyed in all four 
IRAC channels in each individual AOR; W3 Main/(OII), 
KR 140/KR 140-N, and AFGL 333 Spitzer regions, re- 



fectively (Fig. [sj Table 
and 2MASS-based 



§. The fuU list of YSOs (Table 
candidates (Table [2^) appear 
Tn the electronic version of this article. 

For the purpose of the present analysis we chose re- 
liability over completeness, and so this sample requires 
S/N> 5 for any IRAC/MIPS photometry used in any 
particular method. When using 2MASS photometry in 
Phase 2 of the color/magnitude classification scheme, 
those sources with the closest distance to our IRAC 
sources (within 2") were chosen as suitable counterparts 
as long as they had a quality flag 'A', 'B', 'C or 'D' in H 
and K-bands, and at least valid photometry ('D') in the 
J-band (if available) . 

In order to remove as many artifacts and false detec- 
tions as possible, we first bandmerged the IRAC channels 
using channel pairs. A candidate source could still be in- 
cluded in the initial list without suitable photometry in 
all four channels if it either appeared in channels 1 and 
2, or channels 3 and 4. This was done as the first step to 
minimize the need for visual inspection of samples con- 
taining tens of thousands of sources, while at the same 
time attempting to minimize the loss of relatively 'cold' 
sources without a suitable detection at shorter wave- 
lengths. In addition, we also imposed our own 'internal' 
cleaning conditions based on the properties provided by 
APEX and SExtractor for those sources alread y pub- 
lished as reliable detections in the catalog of Ruch et ah] 
( |2007[ ) (e.g., successful PRF fitting, of fitting, clfiptic- 
ity and successful deblending). We note, however, that 
these internal cleaning parameters are physically mean- 
ingless and were only used for 'relative' classification of 
sources within a particular sample, with the only purpose 
being to reject as many artifacts (e.g., PSF residuals and 
spikes around bright sources) and false detections as pos- 
sible. Despite this procedure, these conditions were still 
conservative, and visual inspection and manual rejection 
was still required and performed in the last stages of 
the catalog production process. All sources successfully 
classified as YSO candidates and satisfying all the clean- 



ing conditions and bandmerging requirements have an 
entry= 1 in the 'flag' column in the source catalog (Ta- 
ble [2]). This defines the 'reliable' subset of the final list 
(Fig. §. 

In an attempt to improve the 'completeness' of the 
sample we also used the detections in IRAC channels 
1 and 2 (with the best sensitivity) as a base for a 
new source list. Photometry was performed on fixed 
centroids at longer wavelengths (including MIPS), and 
all the cleaning/selection conditions and the Gutermuth 
classification scheme were again applied to each source. 
Those additional detections satisfying all the catalog re- 
quirements were included in the final catalo g with a flag 
entry= 0. We note that in the scheme of [Gutermuth 
et ai.l (12009'), MIPS photometry was used mainly to re 



classify as YSOs those sources initially rejected as galax- 
ies, AGN, or fake excess sources in the first steps of the 
scheme. 

Photometry derived from the MIPS 24 iim maps was 
generally less reliable. We therefore produced a Cata- 
log 2 (Fig. l6]) based on the same procedure as Cata- 
log 1, but allowing for re-classification of IRAC/2MASS 
sources only if there was a successful MIPS counterpart 
in our original (independently obtained) MIPS source list 
(i.e., the standard fixed-centroid photometry is not used 
in MIPS phase 3). Catalog 2 is therefore more conser- 
vative, because although the SExtractor extraction was 
visually observed to detect all significant sources, a large 
fraction of these detections did not satisfy the cleaning 
conditions after performing APEX PRF fitting on the 
MIPS mosaics, on account of the variable and compli- 
cated background at longer wavelengths. Both catalogs 
(Catalog 1 and Catalog 2) yield very similar source lists 
for Class O/I and Class II sources, differing mainly on 
the number of Class O/I* (highly embedded YSOs) and 
Class II* transition objects. Defined in this work as the 
(*) population, both the highly embedded and transi- 
tion objects rely on MIPS photometry for identification 
and classification. A summary of the number of candi- 
dates found of each class in each field is given in Table 
|4j As expected, highly embedded and transition objects 
are particularly abundant in Catalog 1 due to the use of 
MIPS photometry based on IRAC centroids, with Class 
O/I* sources forming up to ~ 65% of the population in 
W3 Main/ (OH) (- 5.5% Class II*). In Catalog 2 these 
types of objects constitute less than ^ 3.5% of the YSOs 
in each field. 

We created a sample resulting from bandmerging just 
channels 3 and 4 and ran it through the cleaning and 
classification procedures described above. This exper- 
iment produced no new sources. This shows that no 
significant sample of 'cold' sources (without detections 
in IRAC channels 1 and 2) should have been missed by 
the initial bandmerging-by-pairs procedure (within the 
limitations of this technique). We therefore conclude 
that the majority of sources potentially missed in our 
(flag — 1) analysis would have come from the samples 
in IRAC channel 1/channel 2 that do not have counter- 
parts at longer wavelengths due to the diffuse emission 
and sensitivity loss (although we note that Catalogs 1 
and 2 for flag = 1 are very similar; Table j4|. 

Unless mentioned otherwise, in the following sections 
we will use Catalog 1 (all flags) as the primary sample in 
our analysis. As we explain in Section [5j this source list 
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TABLE 2 

YSOS IN EACH SUBREGION OF W3: SAMPLE LIST" 



RA 


Dec 


Class 


Class 


Flag" 


Flag 


h m s (J2000) 


° ' " (J2000) 


Catalog 1 


Catalog 2= 


Catalog 1 


Catalog 2 


34 29 06.92 


61 33 44.75 


classll* 


nomatch 





-1 


34 29 56.38 


61 22 31.69 


classll 


classll 


1 


1 


34 32 05.86 


61 25 22.47 


classll* 


nomatch 





-1 



"Catalog is published in its entirety in the electronic version of this article. 
' Reliability flag (see text): 

1: Candidates satisfying cleaning/bandmerging requirements and with individual 
detections in each band. 

0: Candidates satisfying cleaning/bandmerging requirements but with IRAC 
Channel 3 and 4, and MIPS 24 fim detection (centroid) and photometry 
based on a successful detection in Channels 1/2. 

Same as Catalog 1, but without MIPS 24 /xm- based re-classification 
unless a successful detection was found in our original MIPS source list; i.e., MIPS 
centroid (and photometry) not based exclusively on an IRAC detection. 



TABLE 3 

Photometry for YSOs in each subregion of W3" 



3.6 fim Error 4.5 fj-m Error 5.8 fim Error 8.0 ^m Error 24 /im Error 

12.08 0.01 
12.46 0.01 
12.98 0.01 

"Photometry (magnitudes) for sample in Table|2| 

Catalog is published in its entirety in the electronic version of this article. 



12.05 0.01 11.93 0.03 

12.06 0.01 11.81 0.03 
12.98 0.01 12.86 0.05 



12.03 0.09 9.34 0.19 
11.29 0.08 7.90 0.36 
12.95 0.12 9.03 0.14 




I , , , I , , , I L__L^ L_J! L L__, L_ 

2000 4000 6000 8000 

X [pixel] 



Fig. 5. — Greyscale Spitzer channel 1 mosaic with Class O/I (red), Class II (green) and PMS (blue) candidates. Sample includes all 
sources (and all flags) from Catalog 1, as well as PMS candidates with no IRAC counterparts (see text). 



is expected to be a more reliable indicator of the YSO 
properties. We find that the use (or omission) of the 
(*) population is mostly relevant in those highly popu- 
lated regions with high extinction, mainly IC 1795 and 
KR 140. We cross-correlated our final YSO list with 
the list of infrared sources associated with the cluster IC 



1795 presented in Roccatagliata et al. (2011 1, and we find 
76 YSOs in our sample which are consistent with being 
cluster members. Of these, 30 YSOs belong to the (*) 
population in Catalog 1 (not classified as YSOs in Cat- 



al og 2), and yet all arc also classified as YSOs (Class II) 
in |Roccatagliata et al. (201 1|. When app lying the color 



classification from Megeath et al. ( 2004 ) we find that 



only 4 out of the 76 sources would not have been classi 
fied as YSOs, which supports our decision and the need 
to keep the (*) population in our analysis. Results based 
on Catalog 2 are only mentioned briefly when required. 

We finally note that while the flag = and IRAC short 
wavelength-based catalogs intend to improve the com- 
pleteness of the final sample (and each source was vi- 
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Fig. 6. — Same as Fig[5]but for Catalog 2, without PMS candidates with no IRAC counterpart. 



sually inspected in channel 1), these detections are still 
tentative and should be treated with caution. 



TABLE 4 
YSOs IN THE Spitzer Survey 

W3 Main/ (OH) 
Catalog 1 Catalog 1 Catalog 2 Catalog 2 



Flag 


all" 


1 


all 


1 


Class O/I 
Class O/I* 
Class II 


39 


32 


39 


32 


405 
138 


7 

85 


7 

132 


7 

85 


Class II* 


34 


3 


3 


3 


HAeBe'' 


1 





2 





T-Tauri'' 


110 


50 


126 


58 


KR 140 


Class O/I 
Class O/I* 


88 
130 


80 
3 


84 
3 


78 
3 


Class II 


271 


252 


271 


252 


Class II* 


215 


3 


3 


3 


HAeBe 
T-Tauri 




94 




56 




96 




57 


AFGL 333 


Class O/I 
Class O/I* 


57 
25 


51 



57 



51 



Class II 


140 


129 


140 


129 


Class II* 


24 


2 


2 


2 


HAeBe 














T-Tauri 


86 


43 


86 


43 



" Flag = & Flag = 1 combined. 
^ With Spitzer counterparts. 



We compar ed Catalog 1 and C atalog 2 with the de- 
tections from |Ruch et al.] ( |2007[) for W3 Main/(OH), 
the region analyzed in their study. These authors de- 
tected a total of 295 sources in the four IRAC channels, 
21 (~ 7%) of which were classified as Class I sources, 
and 94 32%) as Class II. All 295 sources were in our 
initial catalog resulting from the SExtractor source de- 
tection process. Our analysis identifies a similar number 
of Class I and Class II candidates in this sample, with a 
total of 117 sources classified as YSOs in Catalog 1 (93 
in Catalog 2). Of the remaining sources in their list not 
classified as YSOs in Catalog 1, 130 are stars, galaxies. 



3 shock/knots of emission or sources with PAH contam- 
inated apertures, and 45 do not satisfy the cleaning and 
reliability conditions. Sources in their sample not classi- 
fied as YSOs in Catalog 2 consist of 149 stars, galaxies, 
3 shock/contaminated aperture sources, and 50 sources 
not satisfying cleaning/reliability conditions. In all cases, 
the percentage of Class II sources is ~ 3 times that of 
Class O/I sources. Table [5] shows the results after apply- 
ing our cleaning conditions and the Gutermuth scheme. 



TABLE 5 

Breakdown of candidat e YSOs with counte rparts in the 
catalog of |ruch et al.| ( |2007| l 



Class 
Flag 



Catalog 
all 



1 



Catalog 1 
1 



Catalog 2 Catalog 2 
all 1 



Total YSO 


117 


77 


93 


77 


ClassO/I 
ClassO/I* 


21 


20 


21 


20 


8 











ClassII 


72 


57 


72 


57 


ClassII* 


16 












Table m includes the number of PMS stars found based 
on 2MASS color and magnitude information, with coun- 
terparts in the Spitzer images. We found 51 more 
T-Tauri stars and 2 HAeBe stars when searching for 
2MASS sources (from the 2MASS PSC) without suit- 
able IRAC counterparts (i.e., IRAC sources satisfying 
our initial cleaning conditions), but located in the com- 
mon (4-channels) areas surveyed by Spitzer. Many of 
these sources are located in the bright, diffuse region sur- 
rounding W3 Main, (OH) and AFGL 333 in the IRAC 
images, which did not allow for proper identification and 
extraction of candidate sources due to confusion. The 
imposed lower limit for the K magnitude of our HAeBe 
stars is within the 2MASS lOcr completeness limit for this 
band (14.3 mag). This sample is also complete in the H 
band ( 15.0 mag), and J band (15.9 mag). The faintest 
2MASS magnitudes for the T-Tauri sample lie beyond 
the 2MASS completeness limits by ~ 1 mag in K and H, 
and ~ 2 mag in J band, and therefore our list will be 
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incomplete at the faint end of the population. 

The use of several steps of cleaning and reliability 
thresholds for our final sample, combined with differ- 
ent sensitivity limits in different regions for a given field, 
complicate the determination of a reliable completeness 
limit for our catalog. In order to provide an estimate 
of the magnitude at which a catalog is 100% complete, 
which due to our cleaning and selection steps differs from 
a fainter limit at which some sources can be detected, we 
examined the number of detections as a function of mag- 
nitude for each field at each IRAC wavelength. Figure^ 
shows the histograms for W3 Main and W3 (OH) for all 
YSOs in Catalog 1, the most complete sample derived in 
this work. Our estimates for 100% completeness, given 
by the 'turnover' points in the YSO distributions for each 
subregion of W3, are included in Table [6) The effects of 
strong, large scale emission in dense regions with high 
stellar activity is evident from the limits derived for the 
HDL. 

Our choice of reliability over completeness excludes 
a large population of faint sources from our final list, 
and therefore the completeness limits for the different 
regions reveal a particularly conservative population for 
IRAC channels 1 and 2. The fiag = sample (using 
the short wavelength channels as a template for pho- 
tometry at longer wavelengths) improves completeness 
at these channels significantly, although we note that our 
estimated completeness limits are still conservative com- 
pare d to samples derived fro m similar Spitzer studies of 
W3 (Polychroni et al. 2010). Clearly, the completeness 
magnitudes quoted tor our work are upper (faint) lim- 
its, especially for the long wavelength channels. The use 
of flag — 1 data or the sample from Catalog 2 would 
be more conservative (and therefore less complete). The 
latter candidate source list, while more reliable than Cat- 
alog 1, is expected to suffer a more severe loss of highly 
embedded sources. This is due particularly to the re- 
jection of poorly fitted sources (failed PRF fitting) with 
APEX, and our strict requirement of suitable fitting and 
photometry in the MIPS band (with severe sensitivity 
loss and confusion due to strong diffuse emission). We 
expect to compensate for this issue with the Herschel 
data currently being analyzed, which will be able to con- 
strain the dust emission much more accurately for those 
embedded (younger) sources. 

With regard to Class types. Class II sources are gener- 
ally weaker at the longer IRAC wavelengths than more 
embedded Class O/I objects, and are therefore more likely 
to be missed b ecause of lower sensitiv ity and bright PAH 
emission (e.g., |Chavarn'a et al.||2008| . 

Despite our attempt to compensate for saturation in 
the IRAC sample by using the short exposure mosaics, 
some sources will still be missed in the brightest and 
most active regions such as W3 Main and W3 (OH) due 
to confusion and location within the bright infrared neb- 
ulosity. The catalogs would also exclude clearly extended 
sources and stellar groups. 

Due to our attempt to provide a reliable sample for 
our PMS sources, we again note that weak-lined and low 
extinction T-Tauri and HAeBe stars have been excluded 
from our catalog in order to avoid major contamination 
from reddened main sequence stars. We cross-checked 
our final list of candidate T-Tauri and HAeBe stars with 
the SIMBAD database, and found no matches with the ex- 



ception of a couple of sources, classified in the database as 
'infrared sources' or 'star in cluster'. Ultimately, spectro- 
scopic surveys of stellar candidates are the most reliable 
method to find and confirm PMS. 

3.2.2. Contamination 

The Spitzer based sample was subjected to strict clean- 
ing conditions and cross-checked with other classification 
schemes to ensure the highest consistency and reliability 
of our sample of YSO candidates. However, it is likely 
that some contamination from extraneous objects like 
galaxies/ AGN, planetary nebulae (PNe) and AGB stars 
is present in the final sample. 

A measure of the contamination by PNe and AGB stars 
was obtaine d by applying the conditions from Robitaille 
et al. ( 2008 ) to the reliable photometry derived from the 
previous YSO classification. For successful classification 
as PN, detections needed to satisfy at least two of the 
four color conditions to ensure reliability in their position 
in the CCDs. While the selection scheme still allows 
for mutual contamination of YSO/AGB stars in both 
samples, this technique can still provide a useful measure 
of the general contamination of our YSO sample. The 
PN contamination using this method is predicted to be 
< 1.5% for W3 Main/(OH), and < 1% and < 2% for 
KR 140 and AFGL 333, respectively. The contamination 
from AGB stars is expected to be < 3.0% and < 0.5% 
for W3 Main/(OH) and AFGL 333. No candidate AGB 
stars are found in KR 140 using this scheme. 

In their analysis, Gutermuth et al. ( 2009j ) adapted their 
classification to work with the larger distances present in 
their survey of up to ^ 1 kpc. The use of this classifica- 
tion at the distance of W3 (~ 2 kpc) is expected to shift 
the proportion of classified YSOs toward the brightest 
sources, because of the loss of dimmer sources rejected 
through the magnitude limits imposed on the sample 
and a possible contamination of the galaxy sample with 
YSOs. To check this effect, we ran the codes on the YSO 
sample in W5 from the work of 'Koenig et al. (2008), ap- 
plying the new classification conditions to the photome- 
try provided by these authors. W5 belongs to a massive 
complex (neighboring W3 and W4) and is considered to 
be at the same distance as these GMCs. These authors 
found a significant proportion of YSOs being misclassi- 
ficd as non YSOs, which was evident, for instance, in the 
'clustered' properties of the 'galaxies'. 
Although the code used in t he present wo r k is a 're- 



vised' version of that used by Koenig et al. (20081, we 



successfully classified 95.5% of their sources as YSOs 
of the same type (excluding the very few sources not sat- 
isfying the requirement of magnitude error < 0.2 mag). 
About two thirds of the 4.5% sources where we disa greed 
with the classification from Koenig et al. (12008) were 
classified as AGN, and 0% as PAM galaxies, nicreasing 
the magnitude limits for AGN classificatio n by 0.5 mag to 



accou nt for the larger distance of W3 (e.g., Megeath et al. 
2009 1 , we find no differences in the percentage of sources 
classified as AGN. Although we cannot perform this test 
on their rejected (non- YSO) sample (as that data was 
not included in their publication), we expect this code to 
detect and recover successfully the main YSO population 
in W3, without the need for major additional modifica- 
tions. 

This is in agreement with the following analysis of 



10 Rivera-Ingraham, A. et al. 



TABLE 6 

IRAC COMPLETENESS LIMITS FOR CATALOG 1 



AOR 


Channel 1 
All/Flag= 1 


Channel 2 Channel 3 
All/Flag= 1 All/Flag= 1 


Channel 4 
All/Flag= 1 


W3 Main/W3(OH) 
KR 140 
AFGL 333 


14.1/12.4 
13.1/13.7 
12.9/12.8 


12.9/11.8 12.7/10.6 
13.1/13.0 13.2/12.8 
12.1/11.8 11.9/11.6 


11.3/10.3 
12.2/11.9 
10.4/10.9 



14 12 10 




16 14 12 10 




Fig. 7. — Number of YSO candidates as a function of magnitude and channel for the region comprising W3 Main and W3 (OH). Channels 
1-4 are shown in order from left to right, top to bottom. Completeness limits for each region (Table[6| have been derived from the turnover 
points of the distributions. 



the distribution of those sources classified as galaxies or 
AGNs. Using the 'Distance to Nearest Neighbor' tech- 
nique ( |Clark fc Evansl[T954| we measured the mean ob- 
served distance between galaxy candidates relative to the 
mean distance that would be expected for a random dis- 
tribution for a population with the same characteristics. 
The ratio between these two quantities (R) would be 
equal to one for a perfectly random distribution. We 
obtain R= 0.94, 0.98, and 0.88, with significance levels 
of - 25%, - 37% and - 5% for W3 Main/(OH), KR 
140 and AFGL 333, respectively. This supports the ran- 
dom spatial distribution of our galaxy candidates and 
low YSO contamination based on the lack of significant 
clustering. Coordinates for the list of galaxy candidates 
is also available online (Table [2]3). 

4. YSO ANALYSIS: OBSERVED AND INTRINSIC 
PROPERTIES 

In this work we aim to produce a reliable list of YSO 
candidates, which will be used in our upcoming papers 
to investigate the early stages of massive star formation 
in W3. An estimate of the mass of each of our YSO can- 
didates in the different classes is a crucial component of 
our analysis. Spectral energy distributions (SEDs) have 
been extensively used to estimate this parameter (e.g., 
Robitaille et al. 2006). However, the highly embedded 
state of some of these sources demands proper modelling 
of the dust envelope and disk, which will be the focus of 
our future Herschel-based analysis. 

Based on t he cla ssification scheme proposed by |Guter^ 
muth et al. ( 2009 1 we have produced a list of YSO can- 
didates in the regions of W3 Main/(OH), KR 140, and 
AFGL 333, without any preliminary bias on selection 
according to clump/core association. This classification 



has been compared to and supplemented by alternative 
classification schemes, as described below. 

We classif ied sources as 'red' acc ording to the nomen- 
clature from Robitaille et al. ( 2008 ), which should include 
all Class I, I' lat, and a la rge number of Class II s ources 
as defined by |Lada| ( [l987[ ) and |Greene et al.| ( [T994[) . Only 
three sources in KR 140 classified as Class 0/1 in Catalog 
1 (and two in Catalog 2) were not initially classified as red 
sources. These sources showed a slight flattening of the 
SED at longer wavelengths which is responsible for this 
misclassification, but all have been visually checked and 
their rising SEDs are consistent with embedded sources. 

In this section we characterize the behavior of the 
Spitzer YSOs in the 2MASS and IRAC color-magnitude 
space. With this analysis we investigate the presence 
of possible identifying characteristics for the different 
classes, as well as near/mid infrared properties that may 
help in the identification of low mass and intermedi- 
ate/high mass YSOs in our sample (e.g., T-Tauri and 
HAeBe PMS stars) based on the information used in this 
work. 

4.1. YSO Stages and the IRAC CCD 



The 'stage' classification from Robitaille et al. (20061 
is particularly useful when combined with the Class 
scheme in order to avoid contradictions between observed 
(color, magnitude) and inferred (e.g., Mcnv, Mdisk) prop- 
erties. Figure |8] shows our YSO sample from Catalog 1 
(flag — 1) and the regions in the CCD that most sources 
of different stages are predicted to occupy (these and all 
the other figures in this analysis do not include error bars 
for clarity). We obtain a similar figure when using the 
candidate list from Catalog 2. While this scheme is useful 
to separate sources with and without circumstellar ma- 
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Fig. 8.— IRAC CCD of ClassO/I+{*) (squares) and ClassII+(*) (diamonds) YSOs from Catalog 1 in W3 Main/(OH ) (left), KR 140 

i middle) and AFGL 333 (right). Black solid lines mark the areas where the majority of Stage I, II and III sources from [Robitaille et al.| 
2006[ l are found. 



TABLE 7 

General 2MASS Properties for 'Class O/I' and 'Class IF Populations 



Class Mean [H-KJ 


(7 


Mean [J-H] 


(T 


Mean K 




Objects 


O/I 


1.5 


0.6 


2.0 


0.8 


13.6 


1.1 


68 


II 0.9 


0.4 


1.4 


0.5 


13.6 


1.1 


412 






T-Test Stat. 


Significance 


F-Test Stat. 


Significance 






O/I vs. T-Tauri 


H-KJ 


11.5 


8E-22 


6.1 


6E-13 






J-H] 


10.0 


3E-18 


9.3 


6E-18 








K] 


-5.0 


2E-6 


1.3 


0.3 






O/I vs. II 


H-K] 


10.8 


2E-24 


2.3 


5E-07 






J-H] 


7.9 


3E-14 


2.6 


lE-08 








K] 


-0.8 


0.4 


1.1 


0.6 






O/I vs. HAeBe 


H-K] 


3.7 


4E-04 


2.8 


4E-03 






J-H] 


6.2 


lE-08 


3.1 


2E-03 








K] 


13.2 


3E-23 


2.2 


0.01 






II vs. T-Tauri 


H-K] 


5.5 


6E-08 


2.6 


2E-06 








J-H] 


7.1 


5E-12 


3.6 


lE-09 








K] 


-5.8 


lE-8 


1.2 


0.5 






II vs. HAeBe 


H-K] 


-2.2 


0.03 


1.2 


0.6 








J-H] 


4.5 


l.OE-05 


1.2 


0.6 








K] 


18.4 


0.00 


2.4 


8E-04 







terial, and a significant proportion of Stage I sources can 
be separated from the remaining population, some Stage 

I objects may however still exist in the regions occupied 
by Stage II (disk domain) and Stage III sources. 

For all catalogs, we observed that more than ~ 95% 
of Class O/I YSOs (flag = 1) and - 90% (flag = 0) are 
classified as Stage I sources. Flag = 1 Class II (including 
Class II*) YSOs are more evenly distributed (~ 50%) 
between Stage I and Stage II, although none of these 
sources is as red in [3.6]-[4.5] as Class O/I (including Class 
O/I*) candidates (Figure J8|. We observe < 1% of Class 

II sources above [3.6]-[4.5[> 1, and none above 1.2. We 
consider this to be the limit separating the area exclusive 
to Class O/I Stage I sources, and the area (bluer colors) 
where Class O/I and Class II Stage I are mixed, perhaps 
indicative of a transition from envelope to optically thick 
disks. 

Very few sources lie in the Stage III area, which con- 
firms the robustness of our sample and ability to sep- 
arate embedded sources and optically thick disks from 
those consistent with photospheric colors and optically 



thin disks. 

4.2. The 2MASS CCD and CMD 

We next analyzed the behaviour of the 'stage' vs 'class' 
classification in the 2MASS CCD and CMD for those 
IRAC sources (Catalog 1, quality flag = 1) with coun- 
terparts in the 21VIASS PSC. As mentioned above, both 
Catalog 1 and 2 yield almost identical samples for Class 
O/I and Class II sources, and therefore the conclusions 
below from this analysis are independent of the catalog 
used. 

4.2.1. Analysis of observed (Class) properties 

In Section [3] we established a scheme to separate the 
populations of T-Tauri and HAeBe stars based on mag- 
nitude and color information. HAeBe stars occupy a dis- 
tinct region in the CCD, but K magnitude data is still 
essential to separate the two samples in the region near 
the T-Tauri locus populated by both types. 

We carried out T-Test and F-Test analyses comparing 
the color and magnitude distributions of various classes. 
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Fig. 9.— 2MASS CCD showing Class O/I IRAC sources with 
2MASS counterparts. Lines and T-Tauri data as in Fig. [s] 



Fig. 12. — Same as Fig. |9] but for Class II sources. 
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Fig. 10.— 2MASS CMD showing Class O/I and Class II IRAC 
sources with 2MASS counterparts. Lines and T-Tauri data as in 
Fig.[4l 



Fig. 13. — Same as Fie. fTol but with HAeBe stars. 
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Fig. 11.— 2MASS CCD showing Stage I IRAC sources with 
2MASS counterparts. Lines and T-Tauri data as in Fig. [3] 

Results are reported in Table [7j Statistically, and as 
in previous cases, the Class O/l, Class II and T-Tauri 
samples are indistinguishable in the 2MASS CMD dia- 
gram, and therefore selection of low mass PMS stars is 
not possible in the infrared without prior knowledge of 
the protostar population. In color, the ClassO/I sample is 
intrinsically redder, with the main population of Class II 
sources lying more intermediate between Class O/I (and 



stage II A 
TTauri )K 




Fig. 14. — Same as Fig. [TT] but for Stage II ; 



closer) to the T-Tauri sample (Figures^ 12 and 10). 
- 90% of Class O/I sources lie between OT <JH-K]<T5, 
compared to Class II sources, located within 0.4 <[H- 
K]< 1.7. While both have similar maximum [H-K] value 
of ^ 3.4, ^ 80% and ^ 50% of Class O/I sources lie above 
[H-K]= 1.0 and 1.5 respectively, compared to ^ 30% and 
~ 9% for the Class II population. 

We find a completely different scenario when focusing 
on the intermediate mass HAeBe population. As shown 
in Figure|4j Herbig stars are not only redder, but brighter 
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than T-Tauri stars, with both populations clearly sepa- 
rated in the CMD. Class O/I candidates are not consis- 
tent with the colors or magnitudes of HAeBe stars (Class 
O/I being redder and dimmer). Class II sources are con- 
sistent in color with Hcrbig stars (i.e., intermediate be- 
tween T-Tauri and Class O/I), but are again dimmer than 
typical HAeBe stars (Figure 13). While a spectroscopic 
analysis is still required to comirm our HAeBe candidates 
as such, we find that, contrary to T-Tauri stars, Herbig 
stars may still be selected without prior knowledge of the 
embedded population. 

4.2.2. Comparative Analysis of Intrinsically Different 
Populations (Stages) 

Using the 2MASS data we find that the populations in 
Stage I and Stage II (as defined in Robitaille et al.|2006 ), 
when treated as a whole, are indistinguish able wit h re- 
spect to observed K magnitude (e.g.. Figs 15 16 and 



171). However, while Stage I sources are founam the 
color region occupied by Stage II, members of the former 
population consistently reach redder colors: ^ 90% of 
Stage I sources have [H-K]> 0.6, - 85% have [H-K]> 0.7, 
- 55% have [H-K]> 1.0, and - 25% have [H-K]> 1.5 
(maximum [H-K]=~ 3.5). The corresponding statistics 
for Stage II are ~ 65%, ~ 40%, 10%, and 2% (maximum 
[H-K]= 2.2). 

F-test significance levels between Stage I-Stage II and 
the T-Tauri population show both groups are statisti- 
cally consistent with having been drawn from the same 
parent population in K magnitude as the T-Tauri sample 
(sig. levels: 0.41 and 0.58, respectively), with Stage II 
also consistent in color (sig. level: 0.39) with the PMS 
population (dominating in the reddening band of typi- 
cal main-sequence and giant stars). Stage I sources show 
a larger scatter in the CCD (e.g.. Figures 11 14), and 



{e.g., Figures 

significance levels obtained from Student's t-tesT~and F- 
test between the two populations agree with both having 
different parent distributions based on color information 
(Table Is]). 

We fmd no clear boundary in the CCD or CMD sepa- 
rating sources with intrinsic, physically different char- 
acteristics, which could easily be explained by differ- 
ent inclination angles. About 90% of Stage O/I sources 
lie in the range 0.5 <[H-K]< 2.0, ~ 90% of Stage II 
between 0.4 <[H-K]< 1.2, leaving the region in [H- 
K]> 1.2 as the area dominated by (younger) sources with 



Menv/M, > 10-^ yr-^ (e.g., Figure 17) 



5. SPATIAL DISTRIBUTION AND CLUSTERING: GROUP 
CLASSIFICATION AND CHARACTERIZATION 

There is strong evidence that star formation in GMCs 
occurs primarily in clusters. The strong link between 
'clustered star- format ion' and massive star formation im- 
plies that in order to investigate the nature and processes 
involved in the formation of the most massive stars it is 
crucial to investigate the characteristics of young stellar 
clusters and their pre-stellar progenitors. This analysis, 
described b elow, was based on our YSO surface density 



maps (e.g., Chavarri'a et al. 2008) and on the so-called 



Minimum Spanning i'ree Algorithm (MST) 



5.1. 



Minimum Spanning Tree Analysis 

In order to identify and characterize regions of YSO 
clustering we first implemented the nearest-neighbor and 




Fig. 15.— 2MASS CMD showing Stage I IRAC sources with 
2MASS counterparts. Lines and T-Tauri data as in Fig. |4] 
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Fig. 16. — Same as Fig. |15[ but for Stage II sources. 
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Fig. 17. — Same as Fig. [15] but for both Stage I and Stage II 
sources. 



MST techniques ( [Gower fc Ross] 1969| . In the MST 
algorithm there is an optimized break length, Dbrcak 
(the branch length whose removal maximizes the num- 
ber of groups, Ng, with members in each group separated 
by distances shorter than this length) for the specified 
NysO: the minimum number of YSOs for the group to 
be considered a 'cluster' (Nyso)- To facilitate compar- 
ison with previous analyses of adjacent regions at the 
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TABLE 8 

General 2MASS Properties: 'Stage I' vs 'Stage IF Population 



Stage 


Mean [H-KJ 


a 


Mean [J-H] 


<j 


Mean K 




Objects 


1 


1.2 


0.6 


1.7 


0.6 


13.6 


1.1 


254 


2 


0.7 


0.3 


1.3 


0.4 


13.6 


1.1 


203 






T-Tost Stat. 


Significance 


F-Test Stat. 


Significance 








[H-KJ 


10.6 


lE-23 


4.1 


4E-23 








[J-H] 


8.0 


lE-14 


2.9 


5E-14 








[K] 


0.4 


0.7 


1.1 


0.7 







same distance (e.g., W5; Koenig et al. 2008) we chose 
Nyso — 10. We also obtained a new set of results for 
Nyso = 5 to investigate the clustering properties at 
smaller scales. Dbroak was estimated from the peak of 
the distribution of the number of groups sati sfying the 
YSO requirement (Ng) as a function of dbroak (Battinelli 
19911; the minimum dbrcak, as well as the incremental 



step, were chosen to be 0.05 pc, similar to the resolution 
of MIPS 24 /im at 2 kpc. 

Dbreak and Ng are expecte d to be affected by incom- 
pleteness and resolution (e.g., Bastian et al.|2007 |. Nev- 
ertheless, this technique is particularly useful when com- 
paring the relative degree of clustering for different types 
of YSOs in different regions with similar data, as in the 
case of the different subregions in W3. Since we can- 
not confirm physical associations in these groups, clus- 
ter membership cannot be determined without additional 
data. Therefore, for the following analysis we define 
the 'sub-branches' resulting for Dbrcak as 'stellar groups', 
to avoid confusion from true clusters and associations. 
Global results for the Spitzer survey are presented in 
Table [9j When a range of break lengths is found to have 
the same maximum Ng we include the range. MST pa- 
rameters for the entire YSO population (analyzed as a 
whole) and for each subregion (each AOR analyzed in- 
dividually) are shown in Tables [9] and 10 respectively. 
Tables also include results for the YSO population as a 
whole and divided into classes. Each Class is also sepa- 
rated according to Nyso hi order to investigate a possible 
hierarchical structure within larger groups. 

Using the group information provided by the MST al- 
gorithm and the technique from [Battinelli 



explored whether there might be a charac 



( 1991 ) we first 



eristic scale in 



the entire Spitzer survey (Tablej9]). 

Our list of YSO candidates (Class O/I, Class IL Class 
O/I* and Class II* candidates: 'AH' class in Table[9]) from 
Catalog 1 yields Dbreak = 0-6 for a minimum group mem- 
bership of Nyso = 10, resulting in 56% of the YSO pop- 
ulation associated with a group. Flag = 1-only sources 
and Catalog 2 (both containing less than half the sources 
in Catalog 1) are more consistent with a larger length and 
grouped fraction of ^ 1.1 pc and 73.5%, respectively. It 
is clear that incompleteness will affect the optimal break 
length, as missing sources will affect the YSO grouping 
and result in spatially larger groups in order to satisfy the 
mini mum group membership requirement. As shown in 
Table [Tol this appears to be an issue mainly when consid- 
ering the (*) population, which constitutes the most un- 
certain sample and the main difference between Catalog 1 
and Catalog 2. Exclusion of this sample affects the Class 
O/I candidate sample in particular, because when the 
highly embedded population is omitted the few remain- 
ing typical Class O/I sources arc more widely distributed, 
resulting in spatially larger groups for a given Nyso • The 



information from Catalog 1 (all flags) is therefore ex- 
pected to be a more (statistically) signiflcant indicator 
of the properties of the overall population, and so this is 
the primary catalog used in the following analysis. We 
also find that all the results derived for Nyso = 5 are 
comparable with the results derived from Catalog 1 and 
Nyso = 10, with a typical Dbroak ~ 0.54 pc containing 
55% of the YSO popul ation. This is in goo d agreement 



with the results from |Koenig et al. ( 2008 ) , who found 



optimal break lengths in the neighboring region of W5 
(for a sample of the same characteristics) of 0.54 pc (YSO 
fraction: 44%). 



5.2. Determination of Group Intrinsic Properties 

Results above indicate that the global YSO popula- 
tion of W3 as a whole shows a tendency to group with 
a scale Dbroak comparable to or less than half a parsec. 
While larger than typical core sizes associate d with the 
formation of individual stars (~ 0.1 pc; e.g., McKee & 



|Ostriker 2007[ Motte et al. 2007), these scale s are con 
sistent with cl ump-like objects 0.5 pc; e.g., Zinnecker 
fc: Yorke||2007 l consider ed to be the likely b irth place of 



stellar clusters. This led Koenig et al. (2008) to conclude 
that these scales might well be typical ot high mass star 
forming regions. However, the relevance of this result and 
its underlying relation to the actual physical processes in 
star formation remains ambiguous. 

It is important to quantify how representative this 
value is of the inter- YSO separations and how relevant 
this grouping is to the original birth configuration and 
conditions of the eventual stellar members. We exam- 
ined the distribution of Dnoarj the distance of a YSO to 
its nearest neighbor, both for the entire sample ( Tabl e 
11) and for members within a specific group (Table 12). 
Typically, D^oar is considerably smaller than half a par- 
sec, and therefore the optimal Dbrcak is more indicative 
of 'inter-cluster' separation than YSO separation, and 
therefore more relevant for cluster formation than that 
of individual stars. We note that embedded clusters of 
massive stars like the one forming IRS5 (W3 Main) or 
the Trapezium in the Orion Nebula have maximum pro- 
j ected stellar separati ons of the order of 0.02 — 0.05 pc 



(Megeath et al. 2005| . These approach the resolution 
limits ot IRAC and MIPS 24 ^m, respectively. There- 
fore, a given Spitzer 'YSO' may in fact contain more 
than one protostar. The clear link between massive stars 
and clusters, however, make the present study a required 
step for understanding the physics behind high-mass star 
formation and the differences with respect to that of low 
mass stars. 

With this goal in mind, we analyzed each subregion of 
W3 in more detail to investigate the underlying proper- 
ties of the stellar groups found by the MST algorithm, 
as well as possible local differences on the intrinsic char- 
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TABLE 9 
Identified Stellar Groups in W3 



All Survey 





Flag= 




Fla 


g= 1 




Catalog 


Class 


Nyso 




Ng 


% Assoc. 


-Dbrcak 


Ng 


%> Assoc. 


1 


AlP 


10 


0.60 


27 


56% 


1.2 


13 


76% 


2 


All 


10 


0.60 - 1.2 


15 


40 - 78% 


1.2 


13 


76% 


1 


ClassO/I 


10 


3.1 - 3.4 


6 


82 - 83% 


3.1 - 3.4 


6 


80 - 82% 


1 


ClassII 


10 


0.85 


14 


55% 


0.85 


14 


56% 


1 


PMS 


10 


0.70 - 1.2 


6 


30 - 61% 








1 


ClassO/I/II 


10 


0.60 - 1.2 


15 


40 - 78% 


0.80 


13 


54% 


1 


All 


5 


0.45 


60 


53% 


0.55 


33 


54% 


2 


All 


5 


0.55 - 0.60 


37 


55 - 60% 


0.95 


33 


54% 


1 


ClassO/I 


5 


2.3 - 2.8 


10 


80 - 88% 


2.3 - 2.8 


10 


80 - 87% 


1 


ClassII 


5 


0.60 


29 


53% 


0.60 


24 


51% 


1 


PMS 


5 


0.65 


14 


37% 








1 


ClassO/I/II 


5 


0.55 - 0.60 


39 


57 - 62% 


0.55 - 0.65 


32 


54 - 63% 



° Includes ClassO/I, ClassO/I*, ClassII, and ClassII*. Other classes exclude 
(*) candidates unless specifically mentioned. 



TABLE 10 

Identified Stellar Groups in Individual Subregions of W3 



W3 Main/(OH) 





Nyso = 


10 




Nyso = 


5 


Class 


-^brcak 


Ng 


% Assoc. 


D break 


Ng 


% Assoc. 


All" 


0.45 


11 


48% 


0.45 


23 


62% 


ClassO/I/II 


0.75 


5 


63% 


0.55 


11 


51% 


ClassO/I 


2.1 - 3.6 


2 


74 - 97% 


2.0 


4 


77% 


ClassO/I+ClassO/I* 


0.45 - 0.55 


6 


35 - 49% 


0.55 


22 


72% 


ClassII 


0.65 - 0.85 


2 


35 - 49% 


0.55 


8 


44% 


ClassII+ClassII* 


0.65 - 0.85 


2 


30 - 43% 


0.55 


9 


40% 


PMS 


0.75 - 1.2 


2 


38 - 85% 


0.65 


8 


45% 



KR 140 


All 


1.1 


14 


67% 


0.85 - 0.95 


30 


72 - 77% 


ClassO/I/II 


0.65 - 


1.6 


7 


42 - 76% 


0.65 


18 


61% 


ClassO/I 


3.1 - 


3.4 


3 


66% 


2.3 - 3.0 


6 


68 - 83% 


ClassO/I+ClassO/I* 


2.2 - 


2.9 


5 


78 - 90% 


0.75 


13 


62% 


ClassII 


0.75 - 


0.85 


8 


54 - 55% 


1.1 


12 


69% 


ClassII+ClassII* 


1.2 


10 


61% 


0.95 


24 


65% 


PMS 


3.9 - 


4.7 


3 


75 - 93% 


2.3 - 2.5 


7 


73 - 77% 



AFGL 333 



All 


0.55 


6 


44% 


0.55 


12 


59% 


ClassO/I/II 


0.55 - 1.3 


4 


41 - 93% 


0.55 


11 


66% 


ClassO/I 


0.25 - 0.35 


2 


40 - 49% 


0.15 - 0.35 


2 


19 - 49% 


ClassO/I+ClassO/I* 


0.25 - 2.7 


2 


28 - 94% 


1.3 - 1.5 


4 


78% 


ClassII 


0.85 - 1.3 


4 


60 - 90% 


0.55 


9 


49% 


ClassII+ClassII* 


0.85 - 1.3 


4 


52 - 85% 


0.55 


10 


46% 


PMS 


0.65 - 1.1 


3 


43 - 60% 


0.65 - 1.1 


6 


60 - 82% 



Class O/I, Class 0/1*, Class II, and Class II* 



acteristics of the steUar population. Tables [TT] and [12] 
include the parameters derived from the YSO candidate 
list in Catalog 1. The former includes, for each subre- 
gion, the total area surveyed, the extinction range, the 
total mass in the region deriv ed from the extinction maps 
(e.g^;^ |Heiderman et al.|r2010| ), the number of YSOs (Ta- 
ble]4 1 , the total surf ace density, the YSO surface density 

and the star formation efli- 



(e.g., Chavarria et al.^^2008) . 



ciency (yi-'Ej ot the region assuming that the YSOs are 
solar-mass stars. The latter assumption is not expected 
to be accurate, especially considering the possibihty of 
some Spitzer YSOs actually being more than one object. 
However, we used this parameter as a measure of the rel- 
ative properties of the different stellar groups, just like 
when considering the 'ages' of the regions in W3 (Section 

t. In both tables uncertainties for the mass and surface 
;nsity have been derived from the statistical uncertain- 
ties in the extinction calculations. These uncertainties 



do not include effects such as variations in the extinc- 
tion law or background fitting uncertainties during the 
creation of the maps. The final errors in these maps 
do not account for the larger differences observed with 
respect to other extinction estimates in the literature, 
and which depen d on knowledge of the du st emissivity 
and temperature ( [Rowles fc Froebrich|2009 ) . Changes in 
these last parameters can, by themselves, affect the esti- 
mated extinction values by a factor of ~ 2, and therefore 
the uncertainties derived from this work will be under- 
estimated. We note that SFE is the amount of mass in 
YSOs at present compared to the total mass, and there- 
fore not nec essarily the ultimate conversion efficiency. 

Table [12] summarizes the average properties of the 
identifiedgroups for Nyso = 10 using the break length 
derived from the global analysis of the entire Spitzer sur- 
vey (Dbreak = 0.6 pc; Table [o]) . Several methods have 
been used to characterize the size and shape of identi- 
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fied groups, from circular to convex-hull tech niques (e.g. 
Bastian et al]|2007[ [Gutermuth et al.||2009[ ). To charac- 



terize the size and elongation ot the group we chose to use 
an elliptical area. The center was defined as the average 
position of the YSOs within the group. The semimajor 
axis is the vector from the chosen center to the farthest 
YSO. The semiminor axis is the minimum size req uired 
to keep all the YSOs within the eUipse (Table [l2]). Pa- 
rameters for individual groups have be en i ncluded in the 
electronic version of this article (Table [l2k) . 
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6. W3 IN PERSPECTIVE: STELLAR CONTENT, CLUSTER 
PROPERTIES, AND STAR FORMATION ACTIVITY 

W3 is be heved to h ave signatures of both triggered 
(HDL; e.K., |Oey e t al.||2005J and isolated (KR140; e.g., 
Kerton et al. 2008") rnassive star formation. The exis- 
tence ot different star formation processes in the same 
cloud makes W3 a prime location for further investiga- 
tion. The identification and characterization of the em- 
bedded/young (1 — 2Myr) cluster population (including 
their relative age) can shed some light on cluster and 
massive star formation. With this goal in mind, in this 
section we use the spatial distribution of YSO classes, 
the properties of the identified groups, and the surface 
density/age maps to present a description of the star 
formation history and activity in W3. Characteristics of 
the groups (including t he (*) population) discussed below 
can be found in Table 12 1 (online). 

We carried out the mam analysis on W3 by working on 
each subregion individually, while using common group 
properties (e.g., Nyso ^-nd Dbroak) for the entire survey. 
Clearly, the definition and properties of the populations 
that are 'clustered' or 'distributed' (objects outside the 
'cluster' boundaries, stars formed in isolation and mem- 
bers displaced from their original birthplaces) depend on 
the fundamental definition of 'cluster' and the chosen 



boundary (see e.g.,|Bres sert et al.|2010 for a compilation 
of recent cluster identification tecnniques) . The applica- 
tion of a common definition and technique throughout 
an entire cloud like W3 (same distance/resolution) can, 
however, be used to determine the relative properties of 
each subregion, avoiding major systematic effects arising 
from distance dependent factors or an assumed defini- 
tion. For W3, this approach is used to investigate the 
possibility of intrinsic differences between stellar groups 
in different environments and with different stellar activ- 
ity (e.g., HDL vs. KR 140), and so illuminate possible in- 
trinsic differences between the clust ered and distributed 
populations (e.g., | Allen et al. 2007). Unless mentioned 



otherwise, the followmg analysis will be based on the 
YSO sample from Catalog 1 and groups with Nyso = 10 
and Dbrcak — 0.6 pc ('red' groups). 

We made use of the YSO density distribution maps 
and carried out an individual analysis of the identified 
groups in each subfield in W3 (W3 Main/(OH), AFGL 
333, and KR 140). 'Age' maps were also created from 
the ratio of the YSO surface density images, for exam- 
ple the ratio of Class II to Class O/I candidates. Peaks 
in this map represent the 'oldest' of those regions con- 
taining YSOs, whereas a low value in a region populated 
by YSOs indicates a relatively young region. The maps 
used in this work are best suited for the analysis of re- 
gions known to contain YSO groups, because the low 
map values in inactive or unpopulated regions depend 
on YSOs rather distant from the pixel in question. Fig- 
ure 18 shows the age map for the entire W3 cloud. The 
'oldest' star forming regions (peaks) are located in the 
KR 140 field, although we note that many of these old 
regions are not always associated with identified groups 
(see below). 

Finally, it is important to distinguish between the ac- 
tual 'age' of the system (from the start of star formation 
activity) and 'apparent age'. In an undisturbed system 
initiated by a short-lived burst of star formation, our 



ratio method yields reliable ages. High density, central 
regions would contain the oldest population and shorter 
free-fall timescales (tg a p~^^'^). Such an idealized sys- 
tem would show a relatively ordered distribution of its 
population, and the age estimate would be an accept- 
able upper limit. On the other hand, our age estimates 
for triggered regions are expected to be biased toward 
younger ages, because there is strong evidence (e.g., IC 
1795; see below) for star formation being an on-going 
process lasting for at least a few Myr, rather than being 
the result of an isolated burst of star formation. Thus, 
on average, the YSO groups in W3 Main/ (OH) might 
only appear to be the youngest using the simple age map 
contours (Fig. 18). 



6.1. Star Formation in the HDL: W3 Main/ (OH) & 
AFGL 333 

AFGL 333 and W3 Main/ (OH) are located in the east- 
ern HDL neighboring W4 (Fig. [l| . The HDL is the most 
dense structure and essentially all the major activity of 
the CMC is found within its boundaries. In both fields 
we have identified signatures indicative of several modes 
of star formation contributing to the overall structure 
and young stellar population. 

6.1.1. The W3 Main/(OH) Field 

The cluster IC 1795 is located at the center of an 
eroded shell-like structure containing the W3 Main and 
W3 (OH) complexes. Both show the most active massive 
star formation, including deeply embedded clusters, H II 
regions, and a trapeziu m-like system riv aling that in the 
Orion nebula (e.g., Meg eath et'ari|2005[) . 

Main/ (OH) 



Figure 19 shows the W3 



i/lOT) field with the 

identified groups, the contours representing the surface 
density distributions for YSOs calculated on a grid iden- 
tical to that of the Spitzer images, and including highly 
embedded and transition stage candidates (Class O/I* 
and Class IF: (*) sample). This figure also includes 'age' 
contours, obtained from the ratio of the surface density 
maps: [Class IRClass H*]/[Class 0/I+Class O/I*]. 

Figure [20] shows the groups found for the same Nyso 
and breaklength as Figure [l9j but excluding the less 
reliable (*) population. For convenience these will be 
specifically referred to as 'magenta' groups, to distinguish 
them from the 'red' groups in Figure |19[ Reduction of 
the number of candidate YSOs yielded smaller and fewer 
groups, as expected. In addition, the omission of highly 
embedded sources, combined with the confusion in the 
regions of strong IR emission, highly affected the num- 
bers of YSOs detected in the innermost and most active 
regions, especially around W3 Main. We find two ma- 
genta groups within red Group 1, concident with IC 1795. 
The oldest (magenta Group 2; Fig. 20) coincides with 



the oldest part of IC 1795. MagentaTJroup coincides 
with a high extinction region and contains numerous H II 
regions and known clusters, including the well known 
maser sources W3 (OH) and W3 (H2O). Although with 
Nm < 25, the characteristics of these two groups are con- 
sistent with those of mid-rich triggered groups in Figure 
19 (see below) linked to massive star formation, with sur- 
face densities greater than 70Mopc~^, Av-mcan > 5.0, 
and mean inter- YSO separations D < 0.25 pc. 
In Figures 21 and 22 we show the groups identified for 
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Fig. 18. — Greyscale Spitzer channel 1 image of W3 with a superposition of 'age' contours from the ratio of Class Il/Class O/I, including 
(*) population. Only specific contours are shown for clarity: relatively old, 3% of map peak value (red); intermediate, 0.5% (magenta); 
and relatively young, 0.05% (green). 

lence are present, to ensure that a gravitationahy (un- 
stable) bound system is formed. 

The privileged location of IC 1795 (Group 1; Fig 
in the most central and dense parts of the HDL, was 
likely key in the formation of its rich population. Al- 
though it depends on the pre-existence of the HDL, its 
location, morphology and population characteristics are 
compatible with isolated/quiescent formation, and not 
necessarily the product of a triggering event as suggested 
in previous studies. A similar (quiescent) origin is also 
suggested for some individual isolated groups in the outer 
boundaries of the field, characterized by a relatively poor 
and yet closely spaced (contrary to the poor groups asso 
ciated with triggered regions discussed above) old popu 



each particular Class (with and without the (*) popula- 
tion) with the relaxed requirement of Nyso = 5- The 
identified groups more closely trace particular structures 
observed in extinction (e.g., filaments), as well as smaller 
separate concentrations of Class O/I and Class II candi- 
dates. 

We suggest that, overall, clustered formation in com- 
pact clumps/cores can be distinguished from the star for- 
mation process associated with young groups that have 
low surface density and extinction and that ar e fo und in 
triggered regions (e.g.. Groups 6 and 7; Fig. 19 1. Fig- 



ure [23 shows the inter- YSO separations as a function of 
SFE. Groups in 'induced' regions with the largest sep- 
arations have been less efficient in forming young stars, 
which could be due to formation in a turbulent environ- 
ment. This could be indicative of the 'distributed' star 
formation mode taking over in those regions that lack the 
conditions, like high column density and low turbulence, 
required to form massive stars and their (richer) parent 
clusters. Indeed, a SIMBAD search reveals no H II regions 
or clusters associated with any of these triggered groups, 
which also lack localized and significant radio continuum 
emission in the Ca nadian Galactic Plane Survey (CGPS; 
Taylor et al']|2003[ ) maps. Rich groups (25 < < 50) 
such as Group U and Group 5 have average inter- YSO 
separations < 0.25 pc and the highest surface densities, 
and both contain within their perimeters a large propor- 
tion of the major massive star activity of the W3 GMC. 
This supports the classical view that one can have clus- 
tered star formation caused by triggering when at early 
stages high surface densities without disruptive turbu- 



lation in a re latively compact configuration (e.g.. Group 



2 in Fig. 19). Group 1 is the next 'oldest' system after 



Group 2, and it is the most massive and richest group, 
with the smallest separation distances down to 0.02 pc. 

W3 Main and (OH) (Fig. [l]) show the highest extinc- 
tion and both are likely to have (at least part) of their 
stellar population ind uced by IC 1795, in agre ement with 



the conclusio ns from Oey et al. (2005) and Polychroni 



et al.| ( |2010[ ). A triggered scenario is supported by i 
an elongated 'older' Class II dominated region within 
Group 1 exten ding in the directions of W3 Main and W3 



(OH) (see Fig. 19); ii) the cavity and shell like structure 



around IC 1795 nesting both systems; and iii) a tendency 
of Class II candidates to be located toward the most cen- 
tral regions, with younger groups of YSOs following the 
shell around Group 1 (in increasing age: Groups 4, 7, 0, 
5, 3, 6; Fig. 19). Using the ratio of Class Il/Class O/I as 
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Fig. 19. — Greyscale Spitzer channel 1 image of W3 Main/{OH) with identified groups for Nyso = 10 ^^'^ Dbroak = 0.6 pc (red ellipses). 
Yellow contours are Class O/I surface density contours between 1 — 5% of peak value ~ 560 YSO pc~^ in 1% steps. Blue contours are Class 
II contours between 5 — 25% of peak value ~ 100 YSO pc~^ in 5% steps. YSO contours have been chosen to span a common YSO range for 
both classes of ^ 5 — 25 YSO pc~^. Green contours are of the ratio of Class Il/Class O/I YSO surface density maps; these include transition 
and highly embedded candidates: * classification. These 'age' contours are for 0.5-2.5% of the peak value of ~ 300 in 0.5% steps, a range 
chosen to highlight the youngest regions. 
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Fig. 20. — Like Fig. |19[ but excluding the less reliable (*) population. Groups are marked as magenta ellipses. Class O/I surface density 
contours (yellow) between 20 — 90% of peak value ~ 8 YSOpc"^ in 10% steps. Class II contours (blue) between 1.5 — 7% of peak value 
~ lOOYSOpc"^ in ~ 1.4%steps. YSO contours have been chosen to span a common YSO range of ~ 1.5 — 7YSOpc~^. 'Age' contours 
(green) are between 1-31% of peak value of ~ 500 in 5% steps. 
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Right Ascension (J2000) 

Fig. 21. — Like Fig. |19| but for Nyso = 5 and D^jcak = 0.55 pc, which is optimal in this r egio n for both Class O/I + Class O/I* groups 
(brown eUipses), and Class II + Class II* groups (blue ellipses). Parameters are from Table [lO] Figure shows the location of the highest 
concentrations of YSOs according to Class. 




Right Ascension (J3000) 



Fig. 22. — Like Fig. |21[ but excluding the (*) population. Now Dbreak = 1.95 pc, whic h is optimal for Class O/I (brown ellipses), and 
Dbrcak = 0.55, which is optimal for Class II (blue ellipses). Parameters are from Table [lO| 



a relative measurement of age, Group 6, the oldest after 
Group 2 and Group 1, happens to be also the nearest to 
IC 1795, while younger groups have larger distances. 

Secondary bursts of star formation, as well as a non- 
negligible star formation duration, are needed to recon- 
cile the characteristics of t he YSO populatio n in IC 1795 
with the age estimated by Oey et al. ( 2005 1(3 — 5 Myr) 
using optical photometry and spectral analysis. We find 



a highly distributed population of Class II and Class 
O/I sources in the vicinity of IC 1795. Star formation 
is estimated to occur in 1-2 dynamical crossing times 



(Elmegreen 20001, but larger age spreads are possible 
it m ulti ple even ts occur within a certain system. Fig- 
ures 21 and 22 show the groups obtained when separat- 
ing the YSOs according to Class with a minimum group 
membership requirement of Nyso = 5. The 'class sub- 
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Fig. 23. — Inter- YSO separations as a function of SFE for W3 
Main/(OH), KR 140, and AFGL 333. YSOs (including highly em- 
bedded and transition candidates) are associated in groups with 
Nyso = 10 s-iid Dbreak = 0.6 pc. Filled symbols mark those groups 
with the youngest ages (Class Il/Class O/I < 1). 

groups' observed within Group 1 and the triggered pop- 
ulation could be indicative of such additional star forma- 
tion events. This combined activity likely reinforced the 
effects of the central cluster when forming the cavity and 
dense surrounding shell hosting the most massive stars 
in W3. 



Contrary to the suggestion of Oey et al. (2005), some 
of the stellar activity in W3 Mam might also have orig- 
inally initiated in quiescent mode. Figure [24] shows an 
'older' (low and intermediate mass) PMS population to- 
ward the outer edge of this region, shown by the ma- 
genta/red contours. On the other hand there are young 
groups at the inner edge of W3 Main, closer to IC 1795. 
I3y calibrating the Class Il/Class O/I ratio of Group 1 to 
the age of IC 1795 we obtain an average age for these 
young groups of ~ 1.5 — 2.5 Myr. In this calculation we 
assumed the onset of star formation occurred in a sin- 
gle event, and so a larger proportion of Class II sources 
implies a more evolved population. We ignored effects 
such as secondary bursts of star formation in the region, 
and assumed that for groups with similar initial stellar 
characteristics there is a direct relation between the class 
ratio and the actual age of the system. A non-negligible 
period of star formation and internal triggering would 
mean that it would be more appropriate to consider our 
age estimates as lower limits. This age estimate for the 
young groups is however comparable to the estimated 
age for the oldest (diffuse) H II regions and the PMS 
population known to domina te in this region (~ lO^yrs; 



Feigelson & Townsley [2008 ). Therefore, triggering (ei 



ther indirectly by IC 1795, and/or by this pre-existing 
PMS population) could have been responsible for the ob- 
served young proto-OB population in the central regions 
of W3 Main, a secondary burst of massive star formation 
in a region with an already highly enhanced surface den- 
sity. Such a scenario would nicely link IC 1795 to one 
of the pos sible models suggested by Feigelson fc Town~ 
ley ( 2008[ ) to explain the origin of W3 Main (option 4 in 
their list), and could explain the anomalous age distri- 
bution (a central young cluster surrounded by an older 
population) described by these authors. 

6.1.2. The AFGL 333 Field 



W3 Main and W3 (OH) dominate both in massive stars 
and in star formation activity. The other HDL field, 
containing AFGL 333, is characterized by a young stel- 
lar population in its most part associated with high ex- 
tinction regions and filamentary structures, that have re- 
mained undetected in the 2MASS-based extinction maps. 
There are also clear similarities and some differences rel- 
ative to W3 Main/ (OH). 

Figures [25] and [2^ show the identified YSO groups 
(with and without the (*) population) for the AFGL 
333 field, including YSO surface density contours and 
age contours. Figures [27] and [28] show the groups identi- 
fied for Nyso = 5 anddifferent classes (each with their 
corresponding Dbrcak)- The oldest and youngest groups 
are clearly separated, with the youngest (brown ellipses) 
localized in the central regions, and the oldest (blue el- 
lipses) toward the outer parts of AFGL 333 (northern 
parts and trigge red regions neighboring W4). 

Group 4 (Fig. [25]) is the oldest and most massive group 
in the field, witnthe largest area, surface density, and 
highest extinction of ~Av-pcak ~ 7.5). From our analy- 
sis we conclude that its population formed first, clearing 
a region and eroding a cavity-like structure much like 
IC 1795 (sizes ~ 6 — 7pc) but in YSOs less populated. 
Within this cavity, the YSO population shows a more 
compact configuration in the eastern side of the group, 
and a more widely separated population toward the west- 
ern side, where some high extinction filamentary struc- 
tures are identifiable. There is also a population of PMS 
candidates toward the edge of this cavity, on the western 
side of Group 4 (opposite W4). Just as with IC 1795, the 
properties of this group (its location in the middle of the 
HDL and the geometrically 'ordered' population, already 
eroding the surrounding material) are compatible with it 
being produced by quiescent formation, rather than hav- 
ing a triggered origin by the Perseus superbubble. 

A significant difference relative to W3 Main/(OH) is 
the abundance of filaments with and without stellar ac- 
tivity. Groups in this field are all relatively poor, with 
only two g rou ps (red Groups 1 and 4) with N> 25 mem- 
bers (Fig. 25) and with the YSO population mainly lo- 
calized to these filamentary structures. The morphology 
of such filaments, seen well in Spitzer channel 4 (Fig. [2]), 
suggests a turbulent origin, a suggestion supported oy 
their remarkable similarity and continuity with nearby 
(non- ext inction) structures. A good example is Group 6 
(Fig. 25), associated with a filamentary structure within 
the cavity caused by the stellar population in Group 4, 
and containing a relatively old YSO population (but still 
younger than the latter). Its associated population, while 
slightly offset from that filament either due to displace- 
ment since formation or clearing of the immediate sur- 
roundings, clearly traces the overall shape of the high 
extinction, indicating a birth association. While the fila- 
ment in Group 6 is easily traced, only dispersed 'rem- 
nants' of high extinction structures are visible at the 
western side of Group 4. This supports the 'older' evo- 
lutionary stage for the latter, where the apparently dis- 
sociated population may be due to the dispersal of their 
common parental filaments. 

A triggered origin is clearly identified for Groups 7 and 
3 (Nm < 25), both in the outer boundary of the HDL 
and associated within structures carved by the activity 
in W4. The pillar associated with Group 7 coincides with 
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Fig. 24. — Same as Fig. |19| with various 'age' contours superimposed. Blue contours are Class Il/Class O/I surface density contours 
between 5-30% of peak value oi ~ 500 in 5% steps, excluding the (*) population. Red and magenta contours are for PMS/[Class 0/I+Class 
II] between 3-10% of peak value of ^ 80 in 1% steps (red), and 10-60% in 10% steps (magenta), with and without the (*) population, 
respectively. Green contours like Fig. |19| 
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Fig. 25. — Greyscale Spitzer channel 1 image of AFGL 333 with identified groups for Nyso = 10 Dtreak = 0.6 pc (red ellipses). 
Yellow contours are Class O/I surface density contours between 5 — 65% of peak value ~ 120YSOpc~^ in 10% steps. Blue contours are 
Class II contours between 3.5 - 43.5% of peak value ~ ITOYSOpc-^ in 10% steps. Crosses are IRAS 02245-1-6115 (bottom) and IRAS 
02244-1-6117 (top). YSO contours have been chosen to span a common YSO range for both classes of ~ 6 — 75 YSOpc"^. Green contours 
are of the ratio of Class Il/Class O/I YSO surface density maps; these include transition and highly embedded candidates: * classification. 
These 'age' contours are for 10-90% of peak value of ~ 300 in 10% steps. Groups like Group 1 are relatively young, while those like Group 
4 are older. 
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Right Ascension (J2000) 



Fig. 27. — Like Fig. [25]but for Nyso = 5, Dbroak = 1-25 pc, and Class O/I + Class O/I* (brown ellipses), and Nyso 
lass ir^fl ^ ' ' ' 



5, Db 



: 0.55, 



and Class II + Class 



(blue ellipses). Parameters are from Table llOl 



a known cluster (IRAS-02252+6120; Bica et al.||2003[ ), 
containing both a Class II population, located toward the 
outer edges of the pillar, and Class O/I sources, mainly 
in the innermost regions. Such a configuration has previ- 
ously been observed in other pillar- like structures (e.g., 
Choudhury et al. [2010 ) and suggested to be the result 
ot radiative driven implosion (RDI) triggered star for- 
mation, in which an ionization/shock front driven by an 
expanding H II region causes a neighboring overdensity 



to collapse, triggering the formation of stars. 

Group 1 is the youngest and the richest group in Fig- 
ure [25] It is located in the innermost regions of AFGL 
333 and it is associated with the most dense and promi- 
nent filament in this field. If filaments are formed (or 
induced to collapse) by turbulence and compression by 
nearby star activity, then Group 1 might have been in- 
duced by the activity at opposite sides associated with 
IRAS 02245-1-6115 and IRAS 02244-^6117, bottom and 
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Fig. 28.— Like Fig. [26] but for Nyso = 5, Dbreak 
(blue ellipses). Parameters are from Table ncn 

top crosses in Figure[25j respe ctively. The latter contains 
bright infrared sources (BIRS; Elmegreen |1980[ ), and the 
former hosts a k nown cluster and massive star activity 
( Bica et al. 2003 ) . These form the brightest regions in the 
mid-intrared m AFGL 333 and both contain a population 
of PMS in their outermost parts. Thus two active older 
regions sandwich the central young stellar population of 
Group 1. 



Right Ascension (J2000) 
: 0.15 pc, and Class O/I (brown ellipses), and Nyso 



5, Dbreak = 0.55, and Class II 



Group 5 (Fig. 25 1, while associated with the same high 
extinction structure as Group 1, is however relatively 
member-poor. Its location away from the influence of 
the main (infrared-bright) star forming activity in AFGL 
333 may explain the low membership, in contrast with 
the richness of Group 1. 

High extinction structures border Group as well (es- 
pecially noticeable in IRAC channel 4; Fig. [2|. The 
associated filaments are located within an evacuated re- 
gion, much like those in Group 6. However, as mentioned 
above, triggering in such regions of low surface density 
will likely result in isolated star formation (instead of 
clustered) . 

The results from the above analysis remain unchanged 
when the (*) population is excluded, except for the omis- 



of column density. We find a similar situation in the 
KR 140 field (below). For small structures and filaments 
prominent in the submillimeter maps, including KR 140- 
N, north of the KR 140 H II region near the center of 
the field, information is lost in the 2MASS-based extinc- 
tion maps. Our Herschel analysis (in preparation) will be 
able to provide accurate masses and properties for each 
of these structures based on temperature-calibrated dust 
emission. 

6.1.3. Modes and Sequential Star Formation 

It has been argued that the presence of massive star 
forming sites along the interface between W3 and W4 is 
evidence that the main activity in the HDL is triggered 
by the expa nsion of the W4 H II region (e.g., , Carpenter] 
et al. 2000 and references therein). This is supported 



by th e estimated ages of the structures (e.g., ]Oey et aL 
20051). We agree that the HDL structure itself was cre- 



261. 



sion of red Groups 2 and 5 (Fig. 

We note that the features seen in extinction in the mid- 
infrared are better tracers of high column density mate- 
rial than the 2MASS-based extinction map, the latter 
missing some m ajor high extinction areas (e.g.. Group I, 
Group 7; Fig. 
sity) estimates 



25] ). Thus the mass (and surface den- 
for associated groups are lower limits. 
Indeed, Group 1 has the largest number of associated 
YSOs, and it is expected to have a large surface density 
comparable to or greater than groups of similar member- 
ship, > 80 M0 pc~^. The extinction structure associated 
with this group is the brightest in the AFGL 333 field 
in the 850 /im SCUBA map, which is an excellent tracer 



ated by W4, with conditions (e.g., turbulence, surface 
density) favorable for star/cluster formation. However, 
whether the formation and ultimate collapse of clumps 
and cores in W3 were directly triggered is less clear. 
Based on the above evidence it seems also plausible to us 
that formation of the HDL was followed in sequence by 
a more quiescent evolution governed by local conditions. 

Local triggering does play a major role in enhancing 
(massive) star/cluster formation: i) internal triggering 
within evacuated regions (inner shells) generates sec- 
ondary bursts of star formation and a distributed popu- 
lation either by forming or collapsing pre-existing small 
overdensities (clumps and filaments; IC 1795, AFGL- 
Group 4, AFGL-Group 0); ii) compression of high density 
regions generates major bursts of star formation, includ- 
ing massive stars and highly embedded clusters (e.g., in- 
ner part of W3-Main, W3 (OH), AFGL-Group 1, AFGL- 
Group 7). 

Overall, the star formation activity and processes in 
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the AFGL 333 and W3 Main/(OH) fields operate simi- 
larly, albeit less vigorously in the former. Environmen- 
tal physical differences between the two regions, such as 
column density distribution and kinematics, will be in- 
vestigated in upcoming papers. The preponderance of 
filaments and the characteristics of these and other star 
forming and starless structures in the HDL will be the 
subject of our Herschel paper currently in preparation. 

6.2. Star Formation in the Central and Western 
Region: KR I4O-N and KR I40 H 11 Region 

The W3 molecular cloud is, overall, in an advanced 
state of evolution. Channel 4 Spitzer images reveal a 
highly turbulent and dynamic environment with a variety 
of structures such as bright rims, pillars with varying 
orientation, cavities, and filaments (Fig. [2]) . In this third 
field, groups have a wide range of ages and properties 
and can be separated mainly into i) groups associated 
with filaments; ii) groups associated with the active star 
forming regions: KR 140-N and KR 140 (Fig[l]). 

Figures I29] and [30] show the YSO groups (with and 
without the (*) population) identified in the central and 
western regions of the W3 GMC. The figures also show 
YSO surface density contours and age contours. Figures 
[31] and [32] are like Figures [29] and [30l but for Nyso = 5 
and separating YSOs accoraing toclass. These figures 
show a tendency for Class II sources to be more widely 
distributed. While this effect could be enhanced by mo- 
tion away from their birthsites, when Class II sources are 
found systematically located preferentially on one partic- 
ular side of the parent structure, like in the case of Group 
2 (Fig. 29 ) , a propagating triggering disturbance seems 
more likely responsible for the observed distribution. 

The Class O/I population is well localized, almost ex- 
clusively in the high extinction filaments and bright rims. 
The Class II population is found distributed through- 
out the entire cloud, following the infrared (PAH) emis- 
sion structures (channel 4) and within the central 'cavity' 
that separates the HDL and the western region compris- 
ing KR 140-N and KR 140. Independent star forming 
events likely produced the observed distributed popula- 
tion. Star formation events throughout the cloud also 
triggered secondary bursts of star formation resulting in 
the eroded and turbulent environment observed at longer 
wavelengths. An almost square lower column density re- 
gion is noticeable in the extinction maps, but we find no 
evidence of any major event that could have caused an 
evacuation. An age sequence could be inferred for the 
star forming regions associated with the HDL. However, 
the structures and groups found i n th e central and west- 
ern regions of the W3 GMC (Fig. 29 ) appear more char- 
acteristic of sequential or even isolated events, perhaps 
more typical of the normal evolution and aging of a cloud 
without a neighbor like W4, and/or perhaps related to 
the original formation of the GMC. 

Groups in the central region of the W3 cloud are 
member-poor (N^ < 25; Groups 6, 5, 2j_10, 3, and 0, 



in order of increasing membership; Fig. 29). Analysis of 
the environment of the filaments (Groups 6, 5, 1, 0, the 



latter also identified as magenta Group 0; Fig. 30 1 sug- 
gests their YSO population could have been triggered by 
external events. The fact that most run parallel to the 
HDL might indicate some interaction (pillars and ele- 
phant trunks in the region point toward the upper (and 



older) part of AFGL 333), but might, on the other hand, 
just reflect the high density there, based on their range 
of ages. Group 10, located in a bright rim of emission 
in the PAH band in the northern parts, does suggest a 
triggered formation. The orientation of other nearby pil- 
lar structures points to a triggering source located in the 
most northern parts of the cloud and outside the Spitzer 
surveyed area. 

The age of the systems in the central regions of W3 
range from the youngest (Group 5 or 10) to the oldest 
(Groups and 3) in the GMC. Group is also the fil- 
ament with lowest extinction. Group 3 exists within a 
cavity-like structure in the neighborhood of weak emis- 
sion pillar-remnants that point toward the location of 
Group 10. This, combined with the lack of association 
with an extinction structure and its low density is sug- 
gestive of an old triggered population that formed from 
material associated with the now low density pillar struc- 
tures. 

The richest groups (N^ > 25; Groups 4, 9, 8, and 2, in 
order of increasing membersh ip; Fig. [29| a,re ass ociated 
with KR 14 0-N, KR 140 (e.g., |Kerton etal.|2001l|Kerton 
et al. 2008), and the major filament east ot this U 11 
region. 

The structure and YSO population of KR 140-N 
(Group 9) support a triggered origin from RDI. The in- 
nermost region of this structure contains a large popu- 
lation of Class O/I sources, with Class II candidates ex- 
tending in front and behind the shocked region. A sim- 
ilar YSO distribution was observed in pillars associated 
with embedded clusters triggered by W4 (e.g.. Group 
7 in the AFGL 333 field; Fig. [25] ) that were also sug- 
gested to have been produced by RDI. The cometary-like 
morphology also resembles that obtained by theoretical 
models of RDI of a cloud ex posed to a high ioni zing flux 
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Bisbas et al. 2011). A sig- 



nificant tail of 'blown' material extends toward the west, 
containing a largely distributed population of Class II 
and some Class O/I sources, the latter generally associ- 
ated with knots of bright PAH emission. Just as Group 7 
in Figure [25] is triggered by W4, the presence of Class II 
sources 'ahead' of the front for KR 140-N suggests trig- 
gering by a source 'external' to this structure and to the 
east. The location of the rim, at the edge of an evacuated 
region, also suggests an external (albeit so far unidenti- 
fied) influence. 

The (low density) YSO population of Group 7 likely 
associated with the shell of KR 140 is the youngest, in- 
dicating active star formation triggered by the ionizing 
star VES 735. We find a population of Class II and PMS 
sources extending toward the north of the H II region 
(Fig. 33 1 that are likely to have originated because of the 
activity in the latter. Group 8, a rich (N,„ > 50) group 
projected on the H II region, contains a mixed popula- 
tion of Class O/I and II sources (ratio ^ 1) indicative 
of an extended period of star formation. This group, to- 
gether with Group 2, are not only the richest in the field, 
but are also associated with the highest surface densities 
(> 60— 7OM0pc~^) and extinctions (Ay-poak ~ 5.5). 
This extinction is in agreement with that estimated for 
the ex citing O star of KR 140, VES 735, by |Kerton etH] 
(19991 from spectroscopy (Ay ~ 5.4) and optical pho- 
tometry (Ay ~ 5.7 ± 0.2), confirming the reliability of 
the extinction map at least within the resolution-related 
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Right Ascension (J2000) 

Fig. 29. — Greyscale Spitzer channel 1 image of the KR 140 field with identified groups for Nyso = 10 a-nd D^rcak = 0.6 pc (red ellipses). 
Yellow contours are Class O/I surface density contours between 10 — 90% of peak value ~ lOOYSOpc"'^ in 10% steps. Blue contours are 
Class II contours between 4.0 — 34.0% of peak value ~ 250YSOpc~'^ in 10% steps. YSO contours have been chosen to span a common 
YSO range for both classes of ^ 10 — 90 YSOpc"^. Green contours are of the ratio of Class Il/Class O/I YSO surface density maps; these 
include transition and highly embedded candidates: * classification. These 'age' contours arc for 4-20% of peak value of ~ 680 in 2% steps. 




Right Ascension (J2000) 

Fig. 30. — Like Fig. |29[ but excluding the less reliable (*) population. Groups are marked as magenta ellipses. Class O/I surface density 
contours between 5 — 95% of peak value ~ 50 YSO pc~^ in 10% steps. Class II contours between 1.5 — 30.5% of peak value ^ 175 YSO pc~^ 
in 10% steps. YSO contours have been chosen to span a common YSO range for both classes of ~ 2.5 — 50 YSOpc^'^. Age contours arc for 
2-52% of peak value of ~ 1000 in 5% steps. 



limitations. 

The filaments associated with Groups 2 and 4 are the 
most prominent of the region. Our analysis indicates 
that Group 2 is in a relatively highly evolved stage and 



was not t riggered by KR 140, in agreement with previous 
analysis ( Kerton et al. |2008 1 . Contrary to the popula- 
tion of Group 4, Class ii sources in Group 2 are displaced 
from the filament toward the north, with a string of Class 



28 



Rivera-Ingraham, A. et al. 




22"^ 20™ 
Right Ascension (J2000) 



Fig. 31. — Like Fig. |29]but for Nyso = 5, Dbroak = 0.75 pc, and Class O/I + Class O/I* (brown ellipses), and Nyso 
and Class II + Class li* (blue ellipses). Parameters are from Table [To] 



5, Db 



: 0.95, 
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Fig. 32.— Like Fig. [30|but for Nyso = 5, Dbroak 
(blue ellipses). Parameters are from Table ncn 



22" 20"" 
Right Ascension (J8000) 

2.25 pc, and Class O/I (brown ellipses), and Nyso 



5, Db 



1.05, and Class II 



O/I still deeply embedded in the innermost regions of the 
filament. The east-west orientation and arc- like shape of 
this structure, as well as the distribution of the Class II 
population and the material traced in the mid-infrared, 
suggest a possible trigger located in the direction of KR 
140-N. Group 4 is younger, with characteristics and ori- 
entation that suggest a link to the activity in KR 140, at 
least in the northernmost parts. 



All member-rich groups are also identified even after 
excluding the (*) sample (Fig. 



30) 



When reducing the 
> 5 we again reach 



membership requirement to Nyso 

similar conclusions (Figs. |31| and [32| 

Overall, Class O/I candidates are confined to the in- 
nermost regions of the filaments (high ellipticity groups), 
and therefore trace these structures with high accuracy. 
We also detect a distribution of mainly Class II groups 
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Right Ascension (J3000) 

Fig. 33. — Greyscale Spitzer channel 1 image of the KR 140 field. Contours like Fig |18[ but with red con tou rs b etween 3 — 33% in 10% 
steps (red). Small blue circles are 2MASS-based PMS candidates. Brown and blue ellipses like in Figs. |29|and|30[ respectively. 



arou nd t he 'cavity-hke' region bounded by Group 10 
(Fig. [29|, KR 140-N, and Group 3. While such a configu- 
ration IS reminiscent of a triggered origin, the distributed 
Class II population might actually be the result of per- 
colating low level spontaneous star formation. A string 
of such Class II groups is observed crossing the entire 
field, extending from the easternmost regions to KR 140, 
and incorpora ting the filaments identified with Groups 



other active (infrared-bright) star forming sites: Group 
10, Group 9 (KR 140-N), and Group 8 (KR 140). Ex- 
cluding the (*) sources we find the youngest group to 
be that associated with the filament closest to KR 140 



2 and 4 (Fig. 31), suggesting that the filaments in the 



southern central and western parts of W3 are the peak 
overdensities in a region of overall enhanced extinction. 

VES 735 might have form ed as part of a growing loose 
association of massive stars. Kerton et al.|(2001) indicate 
that the nearby sources IRAS 02171-h6(teS^nd IRAS 
02174-1-6052 have luminosities consistent with lower mass 
embedded B stars. We find a Class O/I YSO matching 
IRAS 02171+6058 and no significant radio continuum 
emission, which supports their hypothesis of an embed- 
ded B-type star. We do not detect a Spitzer counterpart 
for the other source. A more detailed analysis is required 
to investigate this possible association for VES 735, in- 
cluding the numerous bright infrared stars (BIRS) in this 
regiorp^ 

In Figure [33] we plot the 'age' contours deduced from 
the ratio of Class Il/Class O/I with (*) sources. Peaks in 
the red contours indicate the oldest age, while green rep- 
resents the youngest. The distribution of ages through- 
out the field is evident. Group 7, associated with the 
shell, is the youngest. The age of VES 735 is esti- 
mated to be ~ 2 Myr, c onsistent wi th an H II region o f 



2000), 



1 - 2 Myr ( [Kerton et a l. 1999; Bal lantyne et al. 

and so this age can be considered an upper limit tor the 
YSO population associated with this group. Group 7 is 
followed in age by Group 5, and those associated with the 

The YSO classification and association with clumps and cores 
will be investigated in our upcoming Herschel paper 



(Group 4 in Fig. 29), followed by the population in KR 
140-N (Group 9), the filament associated with Group 2, 
and Group 3 (central W3), with the oldest groups being 
Group 8, projected on KR 140, and Group 0. 

The regions undergoing star formation (Class O/I) in 
the central and western part of W3 are at very local- 
ized spots in the cloud. Neighboring groups have an 
age spread suggestive of individual star formation events 
and evolution. Filaments host some of the oldest pop- 
ulations, while those with triggered morphology (Group 
10, KR 140-N: Group 9) are younger, as expected if actu- 
ally triggered by (and therefore dependent on) a previous 
event. Star formation events must also continue over at 
least a few Myr, as suggested by the Class II population 
distributed throughout the cloud surrounding the active 
star for ming sites. 

Figure 33 also shows the candidate PMS sources, which 
dominate near the KR 140 H II region. Finding the old- 
est population to be associated with KR 140 itself (and 
therefore independent of any previous star formation), 
supports the 'spontaneous' origin of its massive exciting 
star. 

7. CONCLUSION 

By means of Spitzer and 2MASS data we have iden- 
tified and analyzed the young stellar population in the 
W3 GMC. These were classified according to the stan- 
dard 'CLASS' nomenclature for YSOs, and compared 
to other classification schemes based on intrinsic stel- 
lar properties, such as envelope accretion rate and disk 
mass. We find distinct regions in the CCD separating 
the intrinsically (and observationally) young population 
from the disk-dominated sample, with an intermediate 



30 



Rivera-Ingraham, A. et al. 



region likely containing edge-on optically thick disks and 
weaker envelope candidates. 

Observationally, the low-mass PMS population cannot 
be identified unambiguously without first identifying the 
protostellar Class O/I and Class II sources. Intermediate 
mass stars could in principle be identified without the 
need for any previous protostar classification, although it 
is likely that some HAeBe sources would then be missed 
through misclassification as Class II and T-Tauri objects. 

The YSO population was divided into spatial groups 
according to the minimum spanning tree algorithm. We 
also created YSO surface density maps and 'age' maps. 
These data were used to investigate the characteristics 
and history of the star formation in W3. 

The distinctive HDL is no doubt influenced by the ex- 
pansion of W4. The high density conditions there favored 
the formation of particularly rich bursts of star forma- 
tion. The HDL contains the main star formation activity 
of W3 and has signatures of both spontaneous and trig- 
gered star formation from both external (e.g., W4) and 
internal events. Whether W4 was just key to creating the 
initial favorable conditions for (massive) star and cluster 
formation, the very high surface density regions, clumps 
and cores which then collapse at a later stage, or whether 
it also triggered the star formation in such clumps and 
cores is a more subtle question. 

Our finding of a relatively older population in the west- 
ern side of W3 Main and AFGL 333, opposite to shells 
formed by compression by IC 1795 and W4, respectively, 
suggest star formation could have started first in quies- 
cent mode throughout the CMC, including these struc- 
tures as well as the cluster IC 1795 itself. Subsequently, 
triggering mechanisms by the intense activity (e.g., IC 
1795) were responsible for compressing an already dense 
environment, greatly enhancing and/or inducing major 
bursts of massive star and cluster formation (e.g., W3 
Main, W3 (OH)). 

The evolved state of W3 is also evident in the cen- 
tral and western parts of this cloud, an overall highly 
turbulent and eroded environment. This is believed to 
have been produced by numerous individual star forma- 
tion events that were responsible for triggering secondary 
episodes of star forming activity. The central regions be- 
tween AFGL 333 and KR 140 are particularly rich in 
filaments, seen in extinction in the mid-infrared, whose 
formation and star forming activity appear to be asso- 
ciated with this highly turbulent environment. Recent 
star formation is confined mainly to these very localized 
regions. 

The overlapping spatial distributions of the YSO class 
populations in the most active areas of W3 indicate on- 
going periods of star formation in regions of massive 
star and cluster formation, as opposed to a single, major 
short-lived event. 

Based on number of YSO types (in groups) and ratios 
of class surface densities (in regions) as relative measure- 
ments of age, we find that the region comprising W3 
Main and W3 (OH) is the youngest, followed by the 
central/west region of W3 and AFGL 333. However, in 
AFGL 333 confusion does not allow us to detect groups 
in the currently most active areas, and therefore the age 
is representative of the YSO population associated with 
filaments. We cannot determine if the activity in AFGL 
333 did start first or if the activity all across the HDL 



started at similar times triggered by W4, but of all the 
ages derived for the groups associated with isolated fila- 
ments we find those in this region to be the oldest. In- 
deed, the characteristics of the YSO population and indi- 
vidual group ages suggest that the activity immediately 
west of AFGL 333 could have induced some secondary 
triggered events in filamentary structures in the central 
regions of the cloud. 

On-going stellar activity in W3 Main, IC 1795, and 
KR 140 results in younger apparent ages for their associ- 
ated groups. Nevertheless, the KR 140 region appears to 
contain the oldest population in the western part of W3. 
This age is less than that of IC 1795, and the absence 
of a suitable nearby triggering mechanism supports the 
spontaneous origin for VES 735 and stars in this possi- 
ble association. We find that the age of KR 140-N must 
be intrinsically young. This structure shows not only a 
morphology consistent with triggering by RDI, but also a 
well defined distribution of Class II sources surrounding 
the Class O/I population, similar to that associated with 
pillar structures. 

The age of IC 1795 has been estimated to be around 
'-^ 3 — 5 Myr. Even when considering a supersonic ve- 
locity (for typical velocity in the ionized medium) of 
c~ 10km~^, IRS 5 would require ^ 2.5 Myr, and ^ 3 Myr 
for interaction with KR 140-N. If IRS5 itself were actu- 
ally triggered by IC 1795, then this makes the former an 
even more unlikely candidate for triggering activity in 
the western regions. If a source of triggering is indeed lo- 
cated in the HDL, then at least the western side of these 
structures must have been undergoing star formation ac- 
tivity well before the onset of the present major activity 
in the HDL. This is supported by findings of a relatively 
older population in the western side of these structures. 

Results presented above reveal that despite having con- 
sidered a range of relatively evolved stages of pre-stellar 
evolution, the W3 CMC is still a prime target for in- 
vestigating different modes of star formation. We clas- 
sified as a grouped population those YSOs not part of 
a previously cataloged cluster, but belonging to a group 
identified though the MST analysis. 

Primordial/quiescent grouped formation is identified 
by old, relatively isolated systems with the closest inter- 
YSO separations and compact configurations (not neces- 
sarily circular if parental material is filamentary) , whose 
richness depends on the primordial surface den sity of the 



clump/core. An example is Group 2 in Figure 19 Trig 



gering may be a more efficient mechanism for creating 
high surface density structures favorable to high multi- 
plicity and/or for inducing the collapse of such structures 
(e.g., by overcoming internal pressure and turbulent sup- 
port). All current major clustered massive star activity 
in W3 is believed to have been 'triggered' to some ex- 
tent. Quiescent clustered formation might nevertheless 
have occurred initially to produce IC 1795, whose pref- 
erential location in the HDL resulted in a rich popula- 
tion that subsequently triggered secondary bursts of star 
formation within an already overdense region. The out- 
come for the CMC would also depend ultimately on the 
physics, effectiveness, and limitations of the processes 
creating the precursors of groups and clusters in the 
parental cloud (e.g., filaments, clumps, and cores). 

Triggered formation has occurred throughout W3. In 
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some cases, unstable structures with enough surface den- 
sity and low turbulence are created (collect/collapse 
model). Elsewhere there is collapse induced (RDI) in 
the neighborhood of the triggering source. Examples 
for triggered grouped formation are Group 7 in Figure 
[29} a YSO population associated with a shell of mate- 
rial compressed by the expanding KR 140 H II region, 
and Group 9 associated with a cometary like structure 
in KR 140-N. Clustered triggered formation can be ob- 
served in pillar structures such as Group 7 facing W4 
(Fig. 25 1. The present data cannot confirm whether 



(e.g., Group 6, Fig. 25). The fact that the Class II pop- 



pre-existmg cold seeds had formed in the cloud prior to 
the triggering mechanisms, and this issue will be revis- 
ited with the available Herschel data (Rivera-Ingraham 
et al., in preparation). 

Regardless of the detailed sequence of processes, our 
work indicates that triggering is the main mechanism 
associated with those structures dense enough to host the 
current main cluster and massive star formation activity 
(i.e., W3 Main). Clusters formed within such structures 
are the richest. IRS5 in particular, unresolved in our 
data, has been suggested to be a Trapezium-like system 
of proto OB stars within an envelope < 0.02 pc in size; 
its protostellar density, of ~ 0.5 — 5 x l O^pc"'^, makes 



it one of the most dense clusters known (Megeath et al 
2005||Rod6n et al.|2008"] ). This system is itself e mbedHe^ 



in a rnassive cluster of low mass and PM S stars ( [Megeath 
|et al .|1996 Feigelson & Townsley 2008 1 within the region 
of highest extinction in the entire W3 GMC (Ay = 9.6± 
0.3). 

We classified as distributed those YSO candidates not 
in groups and not presently associated with extinction 
structures in the infrared. Some of the distributed popu- 
lation has an origin associated with high extinction struc- 
tures (e.g., filaments, pillars) cleared by the stellar activ- 
ity (e.g., Fig. [251 western part of Group 4). A filament 
can follow a highly irregular and curved morphology, 
likely linked to large scale turbulence in the environment, 
and the YSOs born in the filament will looked highly dis- 
tributed after the disappearance of their parent structure 



ulation of Group 4 is surrounded by possible 'remnants' 
of parent structures suggests lifetimes for the filaments 
and clumps of the order of a few Myr. Whether the star 
formation processes within filaments differ from those in 
clustered formation needs to be determined. Note that 
a population of filamentary origin can be confused with 
that truly originated in isolation. The latter can be ob- 
served in the 'tail' of KR 140-N and in highly turbulent, 
low surface density regions in the vicinity of active star 
formation. An example can be found in the western side 
of W3, where 'knots' of infrared emission trace isolated 
overdensities (clumps, cores), more relevant for the for- 
mation of individual objects. 

Herschel submillimeter imaging data (Rivera- 
Ingraham et al., in preparation) will be used to 
identify and characterize the population of cores and 
clumps in the W3 GMC. By examining how these differ 
in regions of quiescent and triggered activity or with and 
without an associated YSO population, this analysis 
will probe the processes involved in the earliest stages of 
massive star/cluster formation and the different modes 
of star formation. 

In the present analysis we observe a link between the 
star formation modes and their environment. Filamen- 
tary structures appear to trace the turbulence in the 
cloud. Since the environment is a key parameter in un- 
derstanding the processes and properties of the stellar 
progenitors, molecular data (Polychroni et al. in prepa- 
ration) will also be used to characterize the dynamics of 
the material associated with the Herschel sources, and 
link their differing properties back to their environment. 
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